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A fundamental question of physics is what ultimately happens to matter as it is heated or com- 
pressed. In the realm of very high temperature and density the fundamental degrees of freedom of 
the strong interaction, quarks and gluons, come into play and a transition from matter consisting 
of confined baryons and mesons to a state with 'liberated' quarks and gluons is expected. The 
study of the possible phases of strongly interacting matter is at the focus of many research activi- 
ties worldwide. In this article we discuss physical aspects of the phase diagram, its relation to the 
evolution of the early universe as well as the inner core of neutron stars. We also summarize recent 
progress in the experimental study of hadronic or quark-gluon matter under extreme conditions 
with ultrarclativistic nucleus-nucleus collisions. 
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FIG. 1 The phase diagram of H2Q (|ChapIinl . l2007t ). Be- 
sides the hquid and gaseous phases a variety of crystalhne 
and amorphous phases occur. Of special importance in the 
context of strongly interacting matter is the critical endpoint 
between the vapor and hquid phase. 



I. INTRODUCTION 

Matter that surrounds us comes in a variety of phases 
which can be transformed into each other by a change of 
external conditions such as temperature, pressure, com- 
position etc. Transitions from one phase to another are 
often accompanied by drastic changes in the physical 
properties of a material, such as its elastic properties, 
light transmission, or electrical conductivity. A good ex- 
ample is water whose phases are (partly) accessible to 
everyday experience. Changes in external pressure and 
temperature result in a rich phase diagram which, be- 
sides the familiar liquid and gaseous phases, features a 
variety of solid (ice) phases in which the H2O molecules 
arrange themselves in spatial lattices of certain symme- 
tries (Fig. [J). 

Twelve of such crystalline (ice) phases are known at 
present. In addition, three amorphous (glass) phases 



have been identified. Famous points in the phase dia- 
gram are the triple point where the solid, liquid, and gas 
phases coexist and the critical endpoint at which there is 
no distinction between the liquid and gas phase. This is 
the endpoint of a line of first-order liquid-gas transitions; 
at this point the transition is of second order. 

Under sufficient heating water and, for that matter any 
other substance, goes over into a new state, a 'plasma', 
consisting of ions and free electrons. This transition 
is mediated by molecular or atomic collisions. It is 
continuous ^ and hence not a phase transition in the 
strict thermodynamic sense. On the other hand, the 



^ Under certain conditions there ma y also be a true pla sma phase 
transition, for recent evidence see l lFortov et aLl. l2007fl . 
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plasma exhibits new collective phenomena such as screen- 
ing a nd 'plasma oscillations' ( Mrowczvnski and Thomal . 

Plasma states can also be induced by high 
compression, where electrons are delocalized from their 
orbitals and form a conducting 'degenerate' quantum 
plasma. In contrast to a hot plasma there exists in 
this case a true phase transition, the 'metal-insulator' 
transition (jGcb hard. ,1997 : Mott., IMS,). Good exam- 
ples are white dwarfs, stars at the end of their evolution 
which are stabilized by the degeneracy pressure of free 
electr ons ( Chandrasekhaii flOSlHShapiro and Teukolskvl 
I1983D . 

One may ask what ultimately happens when matter is 
heated and compressed. This is not a purely academic 
question but is of relevance for the early stages of the 
universe as we go backwards in time in the cosmic evolu- 
tion. Also, the properties of dense matter are important 
for our understanding of the composition and properties 
of the inner core of neutron stars, the densest cosmic 
objects. Here, the main players are no longer forces of 
electromagnetic origin but the strong interaction, which 
is responsible for the binding of protons and neutrons 
into nuclei and of quarks and gluons into hadrons. In the 
Standard Model of particle physics the strong interaction 
is described in the framework of a relativistic quantum 
field theory called Quantum Chromodynamics (QCD), 
where point-like quarks and gluons are the elementary 
constituents. 

The question of the fate of matter at very high tem- 
perature was first addressed by Hagedorn in a seminal 
paper in 1 965 jHagedornl 1965 ) and later elaborated by 
Frautschi (jFrautschil . 1971 ). The analysis was based on 
the (pre-QCD) 'bootstrap model' in which strongly in- 
teracting particles (hadrons) were viewed as composite 
'resonances' of lighter hadrons. A natural consequence 
of this model is the exponential growth in the density of 
mass states 
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of the constant k is immaterial. The energy that is sup- 
plied is used entirely for the production of new particles. 
This is of course at variance with our present understand- 
ing of the big bang in which the temperature is set by 
the Planck scale T - Mpianck = x/^c/Gn = 1-22 x 10^^ 
GcV where the 'Planck mass' is the mass for which the 
Schwarzschild radius is equal to the Compton length di- 
vided by TT. The quantity Gjv is the Newtonian gravita- 
tional constant and c is the speed of light^ . Referring to 
the Hagedorn paper and th e Friedman model of cosmol - 
ogy, Huang and Weinberg ( Huang and Weinberd . ll97"ol ) 
speculated in 1970 about a limiting temperature also in 
the big bang but noted: Our present theoretical appara- 
tus is really inadequate to deal with much earlier times, 
say when T > 100 MeV. 

The situation changed in the early and mid 1970's af- 
ter it became clear that hadrons arc built from quarks 
and gluons and hence have substructure. In this con- 
text Itoh proposed in 1970 that there might exist stars 
that are entirely made of ve r y ma ssive quarks, rather 
than ordinary baryons ( Itohl . I197G|)'^. The paradox of 
Hagedorn was taken up in 1975 ( Cabibbo and Par"isil 
119751 : IColfins and Perrvl Il97,^ when it was noted that 
the quark-gluon substructure of hadrons opened the pos- 
sibility for a phase transition to a new state of deconfined 
quark-gluon matter, called the 'quark-gluo n plasma'^. In 
close analogy to Fisher's droplet model ( Fished . Il967f ) 
for phase transitions, Cabibbo and Parisi sketched a 
very simple (second-order) phase boundary for the quark- 
hadron transition. They argued that, when matter is suf- 
ficiently heated or compressed, finite-size hadrons begin 
to overlap and quarks and gluons can travel freely over 
large space-time distances. Within this picture, the limit- 
ing temperature Th is in reality close to or even coincides 
with the critical temperature for the phase transition be- 
tween hadrons and quarks and gluons. With point-like 
quarks and gluons the temperature in the early universe 
can grow beyond bounds (big bang singularity). 



This is well verified by summing u p the hadronic st ates 
listed by the Particle Data Group (|Yao et all . \200^ . A 
fit to the data yields Th ~ 160 - 180 MeV. It is then 
easy to see that the logarithm of the partition function 
of such a 'resonance gas' 
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In Z«« (T) = J2 In / dM, 



(2) 

and, hence, all thermodynamic quantities diverge when 
T = Th, which implies that matter cannot be heated be- 
yond this limiting 'Hagedorn temperature'. Here, InZi 
is the logarithm of the partition function for all well iso- 
lated particles with mass m.^. Above a certain mass toq 
all particles start to overlap and from that point on the 
sum is converted into an integral over the mass density 
p(m) and all particles can be treated in Boltzmann ap- 
proximation. For the present argument the explicit value 



II. STRONGLY INTERACTING MATTER UNDER 
EXTREME CONDITIONS 

A. Quantum Chromodynamics 

To understand the salient features of the quark-hadron 
transition and to appreciate the historical developments 
in its physical understanding we need to recall some 
basic facts about the strong interaction. Its mod- 
ern t heory is Quantum Ch romodynamics, introduced in 
1973 (jFritzsch all . \l97^ . This relativistic field theory 
is formulated in close analogy to Quantum Electrody- 



^ In all formulas below we use H = c = 1 

^ At that time quarks were considered very heavy to account for 

the fact that no free quarks were observed. 

* This term was coined by Edward Shuryak l lShurvakl . Il978al ) 
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namics (QED) as a gauge theory of massive fermionic 
matter fields interacting with massless bosonic gauge 
fields. In QED the Lagrangian density for the interac- 
tion of electrons with photons is given by 

Cqed = -^Ff^^Ff''' + ^n^i[d^ + leA^)^ - m^i^ (3) 

where F^^ denotes the field strength tensor of the elec- 
tromagnetic field, which in terms of the vector potential 
Af^ is obtained as 



Fny dnAy — djyAn 



(4) 



The electrons are represented by the four-component 
Dirac spinor field tp of mass me and the electric charge 
e denotes the fundamental coupling constant. The La- 
grangian is invariant under simultaneous (local) gauge 
transformations of the fermion field of the electron and 
the vector potential 



^ ^ e-^x^, A. 



A,. 
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where x{^) is a space-time dependent real valued func- 
tion. The phase factor e~'-^ is an element of the unitary 
group [/(I), which is hence called the 'gauge group' of 
QED. 

Because of the smallness of the fine structure constant 
a = e^/47r ^ 1/137 the evaluation of physical processes 
can be carried out in perturbation theory with high ac- 
curacy (the so calculated value for the magnetic moment 
of the electron agrees, e.g., to experiment within ten dec- 
imals!). Historically this was one of the great triumphs 
of relativistic field theories. 

In QCD, quarks and gluons are the elementary de- 
grees of freedom. Aside from the relativistic quantum 
numbers dictated by Lorentz invariance, quarks come in 
six 'flavors' (up, down, strange, charm, bottom, top). 
To obtain the correct quantum statistics for hadronic 
wave functions it turns out that quarks as well as glu- 
ons also have to carry 'co l or' as an addi tional quantum 
number ( Greenber3 . Il964 iNambiil . Il960f ). The resulting 
Lagrangian density is then given by 

Cqcd = -jG;i,Gr + 9V*(5^+»5sYA;i)g-m,gg (6) 

where q includes the flavor and color quantum numbers 
to be appropriately summed over. The 'strong coupling 
constant' gs is the analog of the electric charge e and 
denotes the quark mass of a given flavor. These masses 
are generated in the electroweak sector of the Standard 
Model via the Higgs mechanism, first introduced in the 
context of superconductivity. We will return to this point 
and its physical implications later. The gauge group 
structure is more complicated than in QED, since three 
colors are required for each quark^. For group theoretical 



consistency also gluons, the force carriers of the strong 
interaction, have to carry color charge (there are eight 
vector potentials A° instead of one). As a physical con- 
sequence they will self- interact. Mathematically this is 
reflected by a modification of the field strength tensor 



(7) 



which now includes a non-linear term. Its form is en- 
tirely dictated by the gauge group (which is now SU (3) 
rather than U{1)) through its 'structure constants' fabc^- 
The group structure is also reflected in the quark-gluon 
coupling through the 'Gell-Mann' matrices Aq which are 
the analog of the SU{2) Pauli matrices. Denoting the 

group of elements SU{i) by J7(x°) = ^ and defin- 
ing A^ = the gauge transformation corresponding 
to Eq. ([5]) now reads 

A^ ^ U{x'')A^U~\x'') + -{d^U{x^))U~\xn- (8) 

It obviously reproduces QED for the gauge group U{1). 

The more elaborate group structure renders QCD 
much more complicated than QED even at the classical 
level of Maxwell's equations.^ 

In any relativistic field theory the vacuum itself be- 
haves, due to quantum fluctuations, like a polarizable 
medium. In QED the photon, although uncharged, 
can create virtual electron-positron pairs, causing partial 
screening of the charge of a test electron. This implies 
that the dielectric constant of the QED vacuum obeys* 
eo > 1- On the other hand, because of Lorentz invari- 
ance, eo/io = 1, i-C. the magnetic permeability iiq is 
smaller than unity. Thus the QED vacuum behaves like 
a diamagnetic medium. In QCD, however, the gluons 
carry color charge as well as spin. In addition to vir- 
tual quark-antiquark pairs, which screen a color charge 
and thus would make the vacuum diamagnetic, the self- 
interaction of gluons can cause a color magnetization of 
the vacuum and make it paramagnetic. This effect ac- 
tually overcomes the diamagnetic contribution from qq 
pairs such that /ig > 1. The situation is somewhat sim- 
ilar to the paramagnetism of an electron gas, where the 
intrinsic spin alignment of electrons overwhelms the dia- 
magnetism of orbital motion. Since /ig > 1 it follows 



^ Quarks form a fundamental representation of the Lie group 

sr/(3). 



® Gauge groups other than ?7(1) were first discussed by Yang and 
Mills in 1954 l|Yang and Milkll95'J ) in the context of 5(7(2) and 
the corresponding field theories are therefore called 'Yang-Mills 
theories'. Since the generators of SU{N) do not commute, such 
theories are also called 'non-Abclian'. 

^ For instance, the wave equation for the vector potentials j4JJ 
is non-linear and its solutions in Euclidean space-time include 
solitons called 'instantons'. 

* Provided the distance r is large enough so that the virtual cloud 
around the test charge is not penetrated. The distance scale is 
set by the inverse Compton wavelength of the electron, which is 
very small. 
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that eg < 1, so that the color-clcctric interaction be- 
tween charged objects becomes stronger as their separa- 
tion grows ('infrared slavery'). In this sense the QCD vac- 
uum is an 'antiscreening' medium. As the distance r — > 0, 
on the other hand, and eg — > 1, and the interaction be- 
comes weaker ('asymptotic freedom'). This gives rise to a 
pronounced variation ('running') of the strong 'fine struc- 
ture constant' = g^/'iir with (space-time) distance or 
momentum transfer Q. Its mathematical form to lead- 
ing order was worked out in 1 973 by Gross and Wilczek 
and independe ntly by Politzer ( Gross and Wilcze3 . ll973l : 
|Politzej . Tl973l) and reads 



127r 



(33-2iV;)ln(QVA|cD) 



<3' > Aqcd 



(9) 

where Aqcd ~ 200 MeV is called the fundamental QCD 
scale parameter. As indicated in Fig. [2] the running of 
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FIG. 2 The running of the fine structure constant of the 
strong interaction with the m omentum tran sfer Q in a col- 
lision of quarks and/or gluons (jBethkd . l2007l ). 

is by now confirmed by experiments to very high preci- 
sion and the authors were awarded the 2004 Nobel Prize 
in physics for their predictions. Even though a mathe- 
matical proof is still missing, it is generally believed that 
the strong increase of the coupling constant for low val- 
ues of Q is responsible for the fact that isolated quarks 
and gluons have not been observed and are permanently 
'confined' in composite hadrons. 



B. Models of the phase diagram 

A simple pictu re of confinement is provided by the 
MIT-Bag model (jChodos et all Il974l ). Here, the ide- 



alized assumption is made that the QCD vacuum is a 
perfect paramagnet with /xg = oo and eg = 0. A hadron 
is formed by carving a spherical cavity (bag) with radius 
R ~ AqJ^j-, « 1 fm out of the physical vacuum. Inside the 
bag the vacuum is trivial, i.e. /ig = eg = 1, and the in- 
teraction between color charges is therefore weak. From 
the boundary conditions on the chromoelectric and chro- 
momagnetic fields it immediately follows that the color 
fields are totally confined within the hadron.^ The cost 
in energy density for creating the cavity is called the bag 
constant B. After filling the bag with three quarks for 
baryons or quark-antiquark pairs for mesons and impos- 
ing appropriate boundary conditions on the quark wave 
functions to prevent leakage of color currents across the 
boundary, B can be determined from a fit to known 
hadron masses. 

For the quark-hadron transition the MIT-Bag model 
provides the following picture: when matter is heated, 
nuclei eventually dissolve into protons and neutrons (nu- 
cleons). At the same time light hadrons (preferentially 
pions) are created thermally, which increasingly fill the 
space between the nucleons. Because of their finite 
spatial extent the pions and other thermally produced 
hadrons begin to overlap with each other and with the 
bags of the original nucleons such that a network of zones 
with quarks, antiquarks and gluons is formed. At a 
certain critical temperature Tc these zones fill the en- 
tire volume in a 'percolation' transition. This new state 
of matter is the quark-gluon plasma (QGP). The vac- 
uum becomes trivial and the elementary constituents are 
weakly interacting since /ig = eg = 1 everywhere. There 
is, however, a fundamental difference to ordinary elec- 
tromagnetic plasmas in which the transition is caused by 
ionization and therefore gradual. Because of confinement 
there can be no liberation of quarks and radiation of glu- 
ons below the critical temperature. Thus a relatively 
sharp transition with 1ST jT^ << 1 is expected. We will 
return to this issue in the section on numerical solutions 
of QCD on a space-time lattice below. A similar picture 
emerges when matter is strongly compressed. In this case 
the nucleons overlap at a critical number density Uc and 
form a cold degenerate QGP consisting mostly of quarks. 
This state may be realized in the inner core of neutron 
stars and its properties will be discussed later. 

In the MIT-Bag model thermodynamic quantities such 
as energy density and pressure can be calculated as 
a function of temperature and quark chemical poten- 
tial^^/ig and the phase transition is inferred via the Gibbs 



The situation is analogous to the case of a cavity in a perfect 
conductor (superconductor) with /i = 0, e = oo except that the 
role of \i and e are interchanged. 

In contrast to water, where the phase diagram is usually char- 
acterized by pressure and temperature, the number density is 
generally not conserved for relativistic systems. Therefore, the 
grand canonical ensemble with state variables temperature and 
quark chemical potential is used. For strong interactions ^lq en- 
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construction of the phase boundary. Under the simpU- 
fying assumption of a free gas of massless quarks, anti- 
quarks and gluons in the QGP at fixed T and one 
obtains the pressure 



PQGp{T,lJ.q) 



37— 
90 



Jfq_ 

27r2 



B 
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To the factor 37 = 16 + 21, 16 gluonic (8 x 2), 12 quark 
(3x2x2) and 12 antiquark degrees of freedom con- 
tribute^^. For quarks an additional factor of 7/8 ac- 
counts for the differences in Bose-Einstein and Fermi- 
Dirac statistics. The temperature dependence of the 
pressure follows a Stefan-Boltzmann law, in analogy to 
the black-body radiation of massless photons. The prop- 
erties of the physical vacuum are taken into account by 
the bag constant B, which is a measure for the energy 
density of the vacuum. By construction, the quark- 
hadron transition in the MIT bag model is of first or- 
der, implying that the phase boundary is obtained by 
the requirement that, at constant chemical potential, the 
pressure of the QGP is equal to that in the hadronic 
phase. For the latter the equation of state (EoS) of 
hadronic matter is needed. Taking for simplicity a gas of 
massless pions of 3 different charge states, which yields 
PTr{T,iiq) = (37r^/90)r'*, a simple phase diagram emerges 
in which the hadronic phase is separated from the QGP 
by a first-order transition line. Taking for the bag con- 
stant the original MIT fit to hadronic masses, B = 57.5 
MeV/fm^ one obtains Tc ^ 100 MeV at = an d 



/ic ^ 300 MeV at vanishing temperature ( Buballal . l2005( ). 

These results imply a number of problems. On the one 
hand, the transition temperature is too small, as we have 
learned. We will come back to this in the next section. 
On the other hand, at 3^^ = fJ-b ^ Mn (mass of the 
nucleon Mjv = 939 MeV), where homogeneous nuclear 
matter consisting of interacting protons and neutrons is 
formed, a cold QGP is energetically almost degenerate 
with normal nuclear matter. Both problems are, how- 
ever, merely of a quantitative nature and can be circum- 
vented by raising the value of B. More serious is the fact 
that, at large fiq, a gas of nucleons because of its color 
neutrality is always energetically preferred to the QGP. 
The biggest problem is, however, that QCD has a num- 
ber of other symmetries besides local gauge symmetry 
which it shares with QED. Most notable in the present 
context is chiral symmetry, which is exact in the limit of 
vanishing quark masses. For massless fermions their spin 
is aligned either parallel (right handed) or antiparallcl 
(left handed) to the momentum. Chirality of a massless 
fermion is a Lorentz-invariant concept, i.e. left(right)- 
handed particles remain left (right )-handed in all refer- 



sures conservation of baryon number and > implies a non- 
vanishing net quark density riq. 

Here it has been assumed that only up and down flavors con- 
tribute significantly to the quark pressure. 



ence frames^^. For physical up and down quark masses 
of only a few MeV this limit is well satisfied when com- 
paring them to typical hadronic mass scales such as the 
mass of the nucleon^'^. Exact chiral symmetry implies 
that only quarks with the same helicity or 'chirality' in- 
teract, i.e. the left-handed and right-handed world com- 
pletely decouple. This means in particular that physical 
states of opposite parity must be degenerate in mass. 

Similar to a ferromagnet, where rotational symmetry is 
spontaneously broken at low temperatures through spin 
alignment, also the chiral symmetry of the strong inter- 
action is spontaneously broken in the QCD vacuum as 
a result of the strong increase of at small momenta 
(Fig. [21). Empirical evidence is the absence of parity 
doublets in the mass spectrum of hadrons. Since mass- 
less quarks flip their helicity at the bag boundary the 
MIT-Bag model massively violates chiral symmetry. For 
the thermodynamic considerations discussed so far this is 
unimportant, but for other aspects of the phase diagram 
chiral symmetry will be crucial. 

There exist effective theories for the strong interaction 
which emphasize the aspects of chiral symmetry and its 
spontaneous breaking in the physical vacuum. One of the 
most thoroughly studied model in connection with the 
phase diagram dates back to early work by Nambu a nd 
Jona-Lasinio (NJL) (jNambu and Jon~ LasinioL [l961allbh 
in 1961, i.e. before QCD was formulated. In its original 
formulation the NJL model was a relativistic field the- 
ory for interacting point-like nucleons of vanishing mass. 
When applied in the context of QCD, the nucleons were 
later replaced by (nearly) massless up and down quarks 
and the model Lagrangian takes the form 



C 



NJL 



^qiird^-mq)q + G[{qqf + iqtj^fq)^] . (11) 



The interaction between quarks and antiquarks is con- 
structed in a manifestly chirally invariant fashion such 
that Cnjl is invariant under left-right transformations 
of the quark fields in the limit niq — s- (chiral limit) 
Gluons do not appear explicitly but are subsumed in an 
effective short-range interaction of strength G between 
the quarks. For sufficiently large G, chiral symmetry 
is dynamically broken in the ground state through the 
condensation of quark-antiquark pairs, i.e. the vacuum 
expectation value (qq) becomes non- vanishing. This is an 
effect that cannot be produced by perturbation theory. 
As a consequence, a gap in the quark energy spectrum 
occurs. This is in direct analogy to metallic supercon- 
ductivity in which, according to the Bardeen-Cooper- 



At the same time massless left- and right-handed fermions trans- 
form into each other under a parity transformation. 
Also the QED Lagrangian JSJ is chirally symmetric in the limit of 
vanishing rrte . On atomic scales this symmetry is however badly 
broken. 

^'^ Since two quark flavors are involved, the transformation group 
is the direct product of the 'isospin group' 5(7(2), acting on left- 
and right handed quarks, i.e. SU{2)i^ X SU{2)j{. 
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Schricffcr (BCS) (iBardeen et all Il957^ theory, pairs of 
electrons interact via the exchange of lattice phonons and 
condense. 

In a quantum field theory the elementary excitations 
of the vacuum are interpreted as particles. In the orig- 
inal NJL model the energy gap determines the mass of 
the nucleon. It is finite even in the absence of a 'bare 
mass'. Thus mass generation becomes intimately linked 
to the non-trivial structure of the vacuum. In particle 
physics this idea of Nambu was new. Replacing nucle- 
ons by quarks, the (nearly massless) quarks acquire a 
'constituent mass' Mg of around 300-400 MeV. Since a 
nucleon consists essentially of three constituent quarks, 
its mass scale is thus explained. It turns out that dy- 
namical mass generation is not only a feature of the 
NJL model but actually happens in QCD itself as can 
be shown from ab-initio solutions of QCD at large cou- 
pling. Figure [3] summarizes the current status of dy- 

J gs com 
masses (Fischer, I2006D 
malism. 



100% 



namical and Hig gs con tributions to the effective quark 
using the Schwingcr-Dyson for- 



0% 




neous rotation of all spins, which costs no energy-"^^. In 
strong interaction physics with two flavors, this mode is 
identified with the pion. The fact that pions are not ex- 
actly massless is related to the finite bare mass of the 
up and down quarks. Nevertheless the pion mass with 
~ 140 MeV is significantly smaller than that of the p or 
the uj meson (~ 800 MeV ~ 2Mq). 

In the 1980's and 1990's the NJL model was used ex- 
tensively in theoretical studies of the phase diagram. 
Since it incorporates spontaneous symmetry breaking 
and the ensuing mass generation, one can address ques- 
tions of chiral symmetry restoration with increasing T 
and fiq and the corresponding medium modifications of 
hadron masses. The quark-antiquark condensate (qq) 
serves as an order parameter for chiral symmetry break- 
ing, analogous to the spontaneous magnetization in a 
spin system. Similar to the Curie- Weiss transition, the 
order parameter vanishes at a critical temperature Tc in 
the chiral limit. This is the point where chiral symmetry 
is restored and the quarks become massless^^. Figure 3] 
displays a prediction for the evolution of the chiral con- 
densate with temperature and quark-chemical potential 
for physical up and down quark masses obtained in mean- 
field theory. 



condensate 
1 



u d 



c b t 



FIG. 3 Fraction of the effective quark mass generated dy- 
namically (light-grey) as compared to that from the Higgs 
mechanism in the electro-weak sector of the Standard Model 
(dark grey). 

As can be seen, the dynamical contribution becomes 
less important the larger the bare or Higgs mass of the 
quark. While the heaviest top-quark mass is entirely gen- 
erated by the Higgs mechanism, for up and down quarks 
close to 99 % of their mass is dynamical. It is thus fair 
to say that almost all of the mass in the visible universe 
is created through the non-perturbative structure of the 
QCD vacuum. 

In QCD, mesons emerge as bound states of quark- 
antiquark pairs with constituent mass. Because of spon- 
taneous chiral symmetry breaking there appears, how- 
ever, a peculiarity that is known from conde nsed matter 
physi cs and was first noted by J. Goldstone ( Goldston^ . 
Il96l[ ). For vanishing (bare) quark mass there must be a 
massless excitation of the vacuum, known as the 'Gold- 
stone mode'. Such highly collective modes occur e.g. in 
spin systems. The ferromagnetic ground state has a spon- 
taneous alignment of all spins. A spin wave of infinite 
wavelength (A ^ oo,fc ^ 0) corresponds to a simulta- 




T [MeV] 



FIG. 4 Evolution of the chiral-condensate ratio {qq}T,fig / (qq) 
plotted on the vertical axis with temperature and quark chem- 
ical potential as predicted by the NJL model (|Heckmannl . 
120071 ). The region of first order transitions, where the con- 
densate ratio jumps discontinuously, is clearly visible. The lo- 
cation of the critical endpoint (CEP) and its projection onto 
the T, pq-plane is also indicated. 



Spin waves obey the dispersion relation E oc k^. In Lorentz- 
invariant theories E on k for massless particles. 



16 



In the NJL model one finds Mg 



■ mg - 2G{qq) 
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While along the T-axis there is a continuous decrease 
indicating a smooth restoration of chiral symmetry, one 
observes along the ^^-axis a first order phase transition 
in which the condensate develops a discontinuity. With 
increasing T this transition becomes weaker and ends in 
a critical endpoint (CEP) where the transition is second 
order. This is analogous to the liquid-gas transition in 
water (Fig. [T]). 



III. RESULTS FROM LATTICE QCD 

As described in the previous section one may predict, 
using the schematic bag model and the NJL model (both 
focusing on different aspects of the strong interaction), 
that upon heating and compression strongly interacting 
matter undergoes a relatively abrupt transition from the 
hadronic phase to the QGP. The relevant scales for this to 
happen are in the realm of very strong coupling as 1. 
Hence, as for the description of any phase transition, the 
application of perturbative methods, which are very suc- 
cessful for QCD processes at high energies, must fail. 
The only known way to solve the QCD equations from 
first principles in the region of strong coupling is to dis- 
cretize the QCD Lagrangian density on a discrete Eu- 
clidean space-time lattice. Here one makes use of the 
formal analogy between Feynman's path-integral formu- 
lation of a quantum field theory in imaginary time t = it 
and the statistical mechanics of a system with tempera- 
ture T = l/r^^. With this method of 'lattice QCD' the 
partition function of the grand canonical ensemble in the 
path integral formulation 



V[A, q] -fo'^ '^''■fv d^^^Qcn-^f^.Q^i) 



(12) 

can be evaluated stochastically via Monte Carlo sampling 
of field configurations, at least at vanishing (see be- 
low). In Eq.in 

the QCD Lagrangian density ([6|) 

From the partition function, the thermodynamic state 
functions such as energy density and pressure can be de- 
termined as 



Cqqj^ denotes the Euclidean version of 



dlnZ 



dT 



N 



d\nZ 



dV 



p = T 

(13) 

in the thermodynamic limit V,N^ oo;N/V = const. 
At least for matter with an equal number of baryons and 
anti-baryons i.e. for vanishing baryo-chcmical potential 



The connection between a quantum system governed by the 
Hamiltonian H and its statistical description is made by consid- 
ering the transition amplitude (/|e~**^|i) from an initial state 
i to the final state /. Comparing this to the partition func- 
tion Z = Tr(e~''^) (/3 = 1/T) one sees that Z can be obtained 
from the transition amplitude by the replacement it = /3, setting 
i = f = n and summing over n. 



fib = '^fJ'q one obtains in this way quantitative predic- 
tions for t he temperatur e dependence of thermodynamic 
quantities (|KarschL[2003 ). Results are displayed in Fig. [5l 
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FIG. 5 The temperature dependence of the energy density 
from numerical solutions of the Q CD equations on a dis- 
crete (Euclidean) space-time lattice (jKarsc hl. l2002h . The red 
and blue lines denote results with two quark flavors (up and 
down) and three flavors (up, down, strange) of equal mass. 
The crosses indicate the realistic case for which the strange 
quark mass is roughly 150 MeV larger than the up and down 
masses. The arrows indicate the corresponding energy den- 
sities and temperatures reached in current (SPS, RHIC) and 
future (LHC) heavy- ion experiments (see sect. 4). For details 
see text. 

To emphasize deviations from the Stefan-Boltzmann 
behavior expected for a free quark-gluon gas one typi- 
cally shows the reduced energy density e/T^ and pres- 
sure p/T^. Near a critical temperature of Tc = 175 MeV 
the reduced energy density shows a rapid {AT/Tc < 0.1) 
variation, which signals the transition from hadronic 
matter to the QGP. The critical energy density e(Tc) 
is 700 ± 300 MeV/fm"^ which is roughly 5 times higher 
than the energy density in the center of a heavy nu- 
cleus like ^°^Pb. At the same time the chiral condensate 
(qq) = dp/dnriq diminishes rapidly near Tc signaling the 
restoration of broken chiral symmetry. As indicated in 
Fig. [5] a systematic discrepancy of about 15 % between 
the calculated energy density (and pressure) and the free 
gas Stefan-Boltzmann limit is observed for T > 2Tc. Al- 
though this is roughly consistent with the first-order cor- 
rection from perturbation theory, the perturbation series 
is poorly conver gent and resu r nmati on techniques have 
to be employed (jBlaizot et al\ . [200^) for a quantitative 
understanding of the high-temperature EoS. 

These ab-initio numerical findings support the simple 
model results for the existence of a QGP transition dis- 
cussed above. In this connection it should be mentioned, 
however, that most lattice calculations still have to use 
unrealistically large values for the light quark masses and 
rather small space-time volumes. With anticipated high- 
performance computers in the range of hundreds of Ter- 
aflop/s, these calculations will be improved in the near 
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future. Ultimately they will also provide definite answers 
concerning the nature of the transition. Among oth- 
ers, this is of importance for primordial nucleosynthesis, 
i.e. the formation of light elements, such as deuterium, 
helium and lithium. In a strongly first-order quark- 
hadron transition, bubbles form due to statistical fluc- 
tuations, leading to significant spatial inhomogeneities. 
These would influence the local proton-to-ncutron ra- 
tios, providin g inhomogeneous initia. 1 conditions for nu- 



cleosynthesis ( Appl egatc e t al\. 19881: Bovanovskv et al\ . 



I2OO6I: ISchwa^ 



1998; Thomas et aL. 119931). Other conse- 



quences would be the generation of magnetic flelds, gravi- 
tational wa ves and the enhanced pro bability of black-hole 
formation ( Bovanovskv et fflII ■ l2006^ . 

At present, indications are that for = 0, relevant 
for the early universe, the transition is a 'cross over', 
i.e. n ot a true phase tra nsition in the thermodynamic 
sense ( Aoki etall . l2006a[ ). Near Tc the state functions 
change smoothly but rapidly, as discussed above. For 
most of the experimental observables to be discussed be- 
low this subtlety is, however, of minor relevance. A cross 
over would wash out large spatial fluctuations and hence 
rule out inhomogeneous cosmic scenarios. Ve ry recent 
studies (|Aoki et all . l2006bl : ICheng et all l2006l ) indicate 
that the exact value of the transition temperature is still 
poorly known. In fact, these investigations have yielded 
values for Tc in the range 150 - 190 MeV. This is in part 
due to difficulties with the necessary extrapolation to the 
thermodynamic (inflnite volume) limit and in part due 
to the general difficulty in providing an absolute scale 
for the lattice calculations. Progress in this area is ex- 
pected with simulations on much larger lattices at the 
next generation computer facilities. 

While at //g = the lattice results are relatively pre- 
cise, the ab-initio evaluation of the phase boundary in 
the (r, /ig)-plane (Fig. [4]) poses major numerical dif- 
ficulties. This is basically related to the Fcrmi-Dirac 
statistics of the quarks and is known in many-body 
physics as the 'fermion-sign problem'. For the integral 
this implies that the integrand becomes an oscilla- 
tory function and, hence, Monte-Carlo sampling methods 
cease to w ork. Only re c ently n ew methods hav e been 
developed (lAUton et all. 120 03: 'Fo dor and Kat3. 120021: 
Ide Forcrand and Philipsen . 2002t iPhilipsenl . 2006t ) to go 
into the region of finite . 

What can be expected? Considering the phase bound- 
ary as a lin e of (nearly) constant energy density , the 
bag model ( Braun-Munzinger and Stachell . Il996a ) pre- 
dicts that the critical temperature decreases with increas- 
ing ^q. By construction the bag model describes a first- 
order phase transition for all chemical potentials. For 
large values of ^q and low temperatures there are indi- 
cations from various QCD-inspired model studies, chiefly 
the NJL model (see Fig. 2]), that the (chiral) phase transi- 
tion is indeed first order. On the other hand, the lattice 
results discussed above seem to indicate that, at very 
small ^q, the transition is a cross over. This would im- 
ply that there is a critical endpoint in the phase dia- 



gram, where the line of first-order transitions ends in a 
second-order transition (as in the liquid-gas transition 
of water). In analogy to the static magnetic suscepti- 
bility xm = dM/dH in a spin system one can define a 
'chiral susceptibility' as the derivative of the in-medium 
chiral condensate (99) t.^ wrt the bare quark mass 
or equivalently as the second derivative of the pressure, 
Xm = d{qq)T,fj.^/dmq = d^p/dm^. Here the quark mass 
niq plays the role of the external magnetic field H . In the 
Curie- Weiss transition xm diverges. The same should 
happen with Xm at the CEP. On the other hand lat- 
tice studies and model calculations indicate that also the 
quark number susceptibility Xn = duq/dfiq = d^p/dfiq 
diverges. This implies that in the vicinity of the CEP the 
matter becomes very easy to compress since the isother- 
mal compressibility is given by kt = Xn/n-q. It is con- 
jectured that the critical behavior of strongly interacting 
matter lies in the same universality class as the liquid- 
gas transition of water (Stcph anov. . ,2004 ). The exper- 
imental identification of a CEP and its location in the 
(T, fiq) plane would be a major milestone in the study 
of the phase diagram. Although very difficult, there are 
sever al theoretical as well as exper imental efforts under- 
way ( Proceedings of Sciencel . [20061 ) to identify signals for 
such a point. For a recent critical discussion concern- 
ing the existence of a CEP in the QCD phase diagram 
see (|Philipsenl . [2007t) . 



IV. EXPERIMENTS WITH HEAVY IONS 



The phase diagram of strongly interacting matter can 
be accessed experimentally in nucleus-nucleus collisions 
at ultrarelativistic energy, i.e. energies per nucleon in 
the center of mass (cm.) frame that significantly ex- 
ceed the rest mass of a nucleon in the colliding nu- 
clei. After first intensive experimental programs at the 
Brookhaven Alternating Gradient Synchrotron (AGS) 
and the CERN Super Proton Synchrotron (SPS), the ef- 
fort is at present concentrated at the Relativistic Heavy- 
Ion Collider (RHIC) at Brookhaven. A new era of ex- 
perimental quark matter research will begin in 2009 with 
the start of the experimental program at the CERN Large 
Hadron Collider (LHC). Here we will not attempt to give 
an overview of the experimental status in this fie l d (for 
recent reviews see (iBraun-Munzinger and Stachell . 120071 : 



iGvulassv and McLerranl . 20051 )1 but concentrate on a few 
areas which in our view have direct bearing on the phase 
diagram. Before doing so we will, however, briefly sketch 
two of the key results from RHIC, which have led to the 
discovery that quark-gluon matter in the vicinity of the 
phase boundary behaves more like an ideal liquid rather 
than a weakly-interacting plasma. 
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A. Opaque fireballs and the ideal liquid scenario 

At RHIC, Au-Au collisions are investigated at cm. 
energies of 200 GeV per nucleon pair. In such collisions 
a hot fireball is created, which subsequently cools 
and expands until it thermally freezes out^^ and free- 
streaming hadrons reach the detector. The spectroscopy 
of these hadrons (and the much rarer photons, electrons 
and muons) allow conclusions about the state of the 
matter inside the fireball, such as its temperature 
and density. The four experiments a t RHIC hav e 

(lAdams et~aZ 



recent ly summarized their results 



2005bl: lAdcox et all 



Back et aU . 120051 ) 



l2005l : 

For a complete 



lArsene et aU . 12005 



overview see 



also the proceedi ngs of the three recent quark matter 
conferences (iProc. Quark-Matter 2005 Conference , 



20061: [pFoc. Quark-Matter 2006 Conferencd . 12007 



Proc. Quark-Matter 2008 Conferencel . l2008^ ■ 




X 



ter follow closely the predictions (iHuovinen et al. I. I2OOII: 
iKolb and Hein j . l2004l : pTeanev et aZ.l . l2002f ) based on the 



laws of ideal relativistic hydrodynamics. By Fourier anal- 
ysis of the distribution in azimuthal angle <& (see Fig. [S]) 
of the momenta of produced particles, the Fourier coef- 
ficient V2 ~ (cos(2$)) can be determined as a function 
of the particles transverse momentum pt- These distri- 
butions can be used to determine the anisotropy of the 
fireball's shape and are compared, in Fig. [7] for various 
particle species, to the predictions from hydrodynami- 
cal calculations. The observed close agreement between 
data and predictions, in particular concerning the mass 
ordering of the flow coefficients, implies that the fire- 
ball flows collectively like a liquid with negligible shear 
viscosity rj. Similar phenomena were also observed in 
ultracold atomic gases of fcrmions in the limit of very 
large scattering lengths, where it was possible, by mea- 
suring 77 through analysis of the damping rates of breath- 
ing modes, to establish that the s ystem is in a strongly 
coupled state (jQ'Hara et ad . l2002l ). 



^s^f^ = 200GeV ^^^Au + 
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FIG. 6 Geometry of the fireball in the plane perpendicular to 
the beam direction in a nucleus-nucleus collision with large 
impact parameter. 

The produced fireball has such a high density and tem- 
perature that apparently all partons (quarks and glu- 
ons) reach equilibrium very rapidly (over a time scale of 
less than 1 fm/c). Initially, the collision zone is highly 
anisotropic with an almond-like shape, at least for colli- 
sions with not too small impact parameter. The situa- 
tion is schematically described in Fig. [5) In this equili- 
brated, anisotropic volume large pressure gradients exist, 
which determine and drive the hydrodynamic evolution 
of the fireball. Indeed, early observations at RHIC con- 
firmed that the data on the flow pattern of the mat- 



A thermal freeze-out is defined as the point in temperature where 
the density of particles with elastic cross section cr becomes small 
enough so that the mean free path A = l/na is larger than the 
system size. 



FIG. 7 The Fourier coefficient V2 for pions, kaons, protons 
and A baryons (with masses of 140, 495, 940 and 1115 MeV, 
respectively) emitted with transverse momentum pt in semi- 
central Au-Au collisi ons at RHIC . The d ata are from the 
STAR collaboration (lAdams et al!. '2005 a|). T he lines cor- 
respond to predictions (|Huovincn et a/.ri200ll ) from hydro- 
dynamical calculations with an equation of state based on 
weakly interacting quarks and gluons. 

This liquid-like fireball is dense enough that even 
quarks and gluons of high momentum (jets) cannot leave 
without strong rcscattering in the medium. This 'jet 
quenching' manifests itself in a strong suppression (by 
about a factor of 5) of hadrons with large momenta 
transverse to the beam axis compared to expectations 
from a superposition of binary nucleon-nucleon collisions. 
The interpretation is that a parton which eventually 
turns into a hadron must suffer a large energy loss while 
traversing the hot and dense collision zone. To make mat- 
ters quantitative one defines the suppression factor Raa 
as the ratio of the number of entries at a given trans- 
verse momentum pt in Au-Au collisions to that in proton- 
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proton collisions, scaled to the Au-Au system by the 
number of binary nucleon-nucleon collisions such that, 
in the absence of parton energy loss, Raa — 1- Corre- 
sponding data are presented in Fig. [H] The strong sup- 
pression observ ed by PHENIX and, i n fact, by all RHIC 
collaborations (lAdams et oil. 2005b Adcox et al\ . 120051 : 
lArsene et all l200,4 iBack et all l2005h demonstrates the 
opaqueness of the fireball even for high momentum par- 
tons, while photons, which do not participate in strong 
interactions, can leave the fir eball unscathed. Theo- 
retica l analysis of these data ( Gvulassv and McLerranl 
I2OO5I : IVitevl . |2006[ ) provides evidence, albeit indirectly, 
for energy densities exceeding 10 GeV/fm'^ in the cen- 
ter of the fireball. Very interestingly, the fireball is 
apparently opaque enough to stro ngly affect th e spec- 
tra of heavy (c a nd b) quarks ( Abelev et all l2007t 
lAdler et aLl . l2006bD . This w a,s not expected in view of 
the a rguments put forward in ( Dokshitzer and Kharzeevl 
l200l[ ). Although the mechanism for heavy quark en- 
ergy loss is not well understood, the data provide evi- 
dence for their scattering and thermalization in the fire- 
ball. This will become important for the disc ussion about 
quarko nia below. There is even evidence ( Adler et all 
l2006al) for the presence o f Mach cone-lik e shoc k waves 
(|Casalderrev-Solana et all l2005l : IStoecked . l200,^ caused 
by supersonic partons traversing the QGP. Apparently 
both elastic parton-parton collisions as well as gluon ra- 
diation contribute to the energy loss but it is fair to 
say that the details of this mechanism are currently not 
wel l understood. The situation is concisely summarized 
in ( Gvulassv and McLerranl . [20051 ) . 




FIG. 8 Resu l ts fro m the PHENIX collaboration (|Adler et all 
l2007l : lAkibal . I2OO6D for the pt dependence of the suppression 
factor Raa- The suppression visible in the data for vr" and 77 
mesons (yellow and red data points) provides evidence for the 
presence of a dense medium scattering partons at high pt and 
degrading their momenta. Photons (purple) which undergo 
only electromagnetic interactions do not exhibit the effect. 
The bands provide an estimate of systematic uncert ainties . 
The y ellow solid line represents a theoretical spectrum (jVitevl . 
I2OO6I ) calculated under the assumption that the initially high 
Pt parton loses energy by gluon radiation in the dense gluon 
gas inside the fireball. 



With the start of the nucleus- nucleus collision program 
at the LHC in the Fall of 2009 the current understanding 
of jet-quenching and of the ideal-fluid behavior of the hot 
fireball will be subjected to decisive tests. At the much 
higher LHC energy, initial temperatures close to 1 GeV 
can be reached and the fireball is probed with partons 
in the 100 GeV range. It will be exciting to see how the 
currently developed concepts will evolve with the data 
from this new era. 



B. Hadro-Chemistry 

In ideal hydrodynamics no entropy is generated dur- 
ing the expansion and cooling of the fireball, i.e. the 
system evolves through the phase diagram along isen- 
tropes, starting in the QGP phase. This can be experi- 
mentally verified through the production of a variety of 
mesons and baryons. The analysis of particle production 
data at AGS, SPS and RHIC ene r gies has clearly d emon- 
strate d dAndronic et all l2006t feecattini etai\ . |2004 
I2006I : Ib raun-Munzinger et all l2004al ) that the measure- 
ments can be understood to a high accuracy by a sta- 
tistical ansatz in which all hadrons are produced from a 
thermally and chemically equilibrated state. This hadro- 
chemical equilibrium is achieved during or very shortly 
after the phase transition and leads to abundances of the 
measured hadron species that can be described by Bose- 
Einstein or Fermi-Dirac distributions 



/•CSC 



(14) 



of an ideal relativistic quantum gas. Here = M'j + 

Pj^ is the relativistic energy- momentum relation of free 
hadrons of mass A/j, /ij the chemical potential of this 
species, and gj counts the number of degrees of freedom, 
such as spin and charge state of a given hadron. The 
results of such an analysis for the measured abundances 
in central Au-Au collisions at a cm. energy per nucleon 
pair of y/SNN = 130 GeV at RHIC are shown in Fig. O 

Such calculations give, for each beam energy, a set 
of two thermodynamic variables, namely temperature T 
and baryo-chemical potential fj,b at the point of hadropro- 
duction, i.e. at chemical freeze-out^^. This is consistent 
with the assumption that all particles were produced at 
the same instant, i.e. at the same temperature and chem- 
ical potential. Such analyses also provide a striking con- 
firmation for the concept of a lim iting temperature Th 
discussed above ( Hagedornl . Il965[ ). as shown in Fig. [TUl 

The significance of these results is further appreciated 
by entering the (T, fib) values of fixed beam energy into 
the phase diagram (Fig. [TT|) . establishing the 'chemical 



Chemical freeze-out occurs, when inelastic collisions between 
particles cease such that the abundance ratios do not change 
anymore. 



11 



10 



10 



10 



=130 GeV 



• Data 

Model 

T=165.5, fi|,=38 MeV 



n K* K 



n (|) K 



7C Jl 71 7C 



K K 



FIG. 9 Ratios of particle multiplicities in central Au- 
Au collisions at RHIC energies in comparison with a 
fit (|Andronic et ali |2006| ) from a statistical model of ther- 
mally and chemically equilibrated hadrons. 



freeze-out curve' (iBraun-Munzinger and Stachell . Il996bl 
119981: iBraun-Munzinger et a^l.ll996^ . It was noted early- 
on (jClevmans and Redlich . 1998f) that this curve can 
be understood phenomenologically by assuming that the 
freeze-out takes place at a constant energy per particle 
of about 1 GeV. 

In Fig. [11] the experimental points for chemical freeze- 
out ar e compared with the ph ase boundary from lattice 
QCD (iFodor and Katd. l2004^ and fro m the bag model 
( Braun-Munzinger and Stachell . [T996al ). For illustration 
also a recent theo retical prediction for th e possible loca- 
tion of the CEP (|Fodor and Katd . |2004) is shown (tri- 
angle) as well as the trajectory that the early universe 
has taken in the standard big bang model.^" It is in- 
teresting to note that, for ^j, < 250 MeV, the experi- 
mental freeze -out points are close to the calculated phase 
boundary of ( Fodor and Katd . 12004 ). This does not come 
as a surprize. On the contrary, there are good argu- 
ments that the phase transition itself is responsible for 
the e quilibration of all hadrqn specie s. In a recent anal- 
ysis ( Braun-Munzinger et al\ . l2004bl ) it was shown that, 
because of the strongly increasing particle density near 
the phase boundary, multi-particle collisions dominantly 
contribute to the particle production. This leads to a 
rapid equilibration (r < 1 fm) of A baryons and even of 
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FIG. 10 Energy dependen ce of thermal paramete rs T and ji^ 
from a statistical analysis (jAndronic et alV l2006t ) of hadrons 
produced in central nucleus-nucleus collisions. 



baryons with multiple strangeness content (S, baryons) 
compared to the typical expansion time of the fireball of 
several fm. It also explains naturally why all particles 
freeze out in a relatively short time at nearly constant 
temperature. Similar conclusions, alt hough based on dif- 
ferent arguments, can be fou nd in ( Heinz and Kestinl . 
120061: iHeind . Il998t |Stoc3, [199^. The general result from 
these findings is that, at least at small /i^, the temper- 
atures extracted from the chemical analysis are closely 
linked to Tc obtained from the calculated QCD phase 
boundary. Thus, for the first time, a fundamental param- 
eter of the phase diagram, namely the critical tempera- 
ture Tc at small /it., has been confronted with experiment 
and the agreement is very good^^. 

For larger values of jib the measured freeze-out points 



This trajectory is evaluated llBraun-Munzinger and Wambacbl . 
l2006tl by assuming that the early universe expanded isentrop- 
ically under the conditions of charge neutrality and net lep- 
ton - net baryon number conservation and that the entropy per 
baryon is fixed using the known baryon to pho t on ra tio, see also 
llFromerth and Rafelskil 120021: iKampfer et ali l2007t) . Further- 
more, the evolution proceeds in full chemical equilibrium among 
hadrons and leptons until the neutrinos freeze out at a time of 
about t = 1 s. 



We neglect here the above disc ussed uncertainty in Tc obtained 
from recent lattice calculations l|Aoki et al ]. l2006bl : [Ch eng et all 
120061) . It should be pointed out, however, that a value of Tc = 
190 MeV is inconsistent with the scenario discussed here and 
would probably imply the presence of an ultra-dense hadronic 
phase between chemical freeze-out and the phase boundary. 
There is currently no indication of such a phase from experi- 
ment. 
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C. Medium modifications of vector mesons 
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FIG. 11 The QCD phase diagram and data from a chemi- 
cal freeze-out analysis in nucleus-nucleus collisions at various 
beam energies. The data points closest to the T-axis are from 
the highest collision energies. For the other entries see text. 



deviate from the predictions of lattice QCD (Fig. [TT|) . 
At present it is hotly debated whether this deviation 
indicates the existence of a highly compressed hadronic 
phase between the QCD phase boundary and the chem- 
ical freeze-out line, or whether the calculation of the 
phase boundary at large /ib will be modified by signif- 
icant corrections from realistic quark masses and larger 
space-time lattices. Important new insight is expected 
from measurements with the 'Compressed Baryonic Mat- 
ter' (CBM) experiment planned at the future 'Facility for 
Antiproton and Ion Research' (FAIR) at GSI in Darm- 
stadt, as well as from improved lattice simulations. 

The thermalization described above implies that equi- 
librated matter is produced in high-energy collisions be- 
tween nuclei. In e+e^ and hadron-hadron collisions, 
such an equilibration, i n particular in the strangeness 
sector, is not observed ( Braun-Munzinger et al\ . l2004al ) 
although thermal f eatures are observed in the yields of 
prod uced particles ( Becattin"il . ll996l : lBecattini and Heinj . 
I1997D . For very recent discussions of differences and 
simi larities between e^e~ and nucleus- nucleus co llisions 
see ([Andronic et al\ . l2008al : iBecattini efo/] . I2008D . 



For the high-energy domain accessible with Pb ions 
at the LHC the scenario described implies essentially 
small changes in hadron production (apart from an over- 
all yield factor due to the much larger volume). Any 
deviation would be a major surprize and would likely 
indi cate new physics. For specu lations in this direction 
see ( Rafelski and Letessierl . [20081 ) . 



As the spontaneously broken chiral symmetry of the 
strong interaction gets restored at high temperatures and 
large chemical potentials, the quarks loose their 'con- 
stituent' mass and only the 'bare' masses generated in 
the Higgs sector of the Standard Model are left. As 
seen from Fig. [3] this effect is most dramatic for up and 
down quarks and to a somewhat lesser extent also for 
strange quarks. Most naively the mass of a hadron is a 
multiple of the constituent quark mass Mq (for baryons 
M\, ^ 3Mq and for mesons Mm ~ 2Mg) and one would 
therefore expect that all hadron masses consisting of light 
u,d and s quark s should decr ease significantly near the 
phase boundary ( Brownl . ll988l ). More general arguments 
along these lines led to the conjecture of a general scaling 
law in which (nearly) all light hadrons consisting of u,d 
quarks change with some power of the chiral condensate 
ratio (Brown an d Rho, 1991) ('Brown-Rho scaling')^^ 



Mh^{(qq)T,^J{qq)r 



(15) 



Another obvious source of medium modifications of 
hadrons is the increased collision rate in a hot and dense 
medium. As a consequence, many new decay channels 
open, resulting in large widths. Finally, based on chi- 
ral symmetry alone and its spontaneous breaking in the 
vacuum, it can be argued that the spectral properties 
of hadrons with opposite parity become more and more 
similar as the chiral phase transition is approached. 

Since possible modifications of hadron properties 
(masses, decay modes) occur in the hot and dense phase 
of a heavy-ion collision, one needs an experimental probe 
that is sensitive to this state of the matter. More than 
30 yea rs ago it was suggested ( Feinber j . Il976l : IShurvakl 
Il978bl ) that real or virtuaP'^ photons arc ideal, since they 
interact only electromagnetically with the surrounding 
matter and hence leave the reaction zone almost undis- 
turbed. Even at the highest temperatures and compres- 
sion reached in relativistic heavy-ion collisions the mean 
free path of photons is typically 10^ — 10"* fm, which is 
much larger than the size of the fireball. 

Both longitudinal and transverse photon polarizations 
contribute to the di-lepton rate, while real photons can 
only be transversely polarized. According to Fermi's 
Golden Rule the production cross section is directly re- 
lated to the (auto)correlation function {jei,njeim) '^^ 
electromagnetic current which involves the charge carri- 
ers of the system. Taking quarks as fundamental con- 
stituents of strongly interacting matter, jeim is given by 



Jelr. 



E 



The pion is special because of its 'Goldstone character' and there- 
fore its mass should remain largely unaffected. 
Virtual time-like photons correspond to the process of di-lcpton 
(e+e~ or fi'^fi~) pair production or annihilation. 
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3 ' 



(16) 



(For the measurements discussed below only the three 
hght quark flavors are relevant). It is well estabhshed 
by precision measurements that the e+e~ annihilation 
cross section below cm. energies of ~ 1.2 GeV is essen- 
tially saturated by the light vector mesons p, w, (f> with 
the p- meson giving the largest contribution (~ 9 : 1:2). 
Therefore, the in-medium modification of the p-meson in 
di-lepton production in heavy-ion collisions is of particu- 
lar interest. Also the large (F = 150 MeV) width implies 
that the /9-meson decays and is regenerated several times 
during the lifetime of the fireball: the resulting di-leptons 
then carry information about its interior. 

In physical terms the di-lepton signal is. therefore, 
dominantly due to pion annihilation tt^tt" ^ p — s- e~^e~ 
in the hadronic phase or quark annihilation qq e^e~ in 
the partonic phase. If we assume that the fireball formed 
in an ultra-relativistic nucleus-nucleus collision is close to 
thermal equilibrium then the above formalism leads to 
di-lepton (photon) spectra after convolution of the rele- 
vant transition rates with the hydrodynamic space-time 
evolution of the system. 

Experiments to measure di-lepton production in nu- 
clear collisions have been conducted since the late 
1980's starting with data taking at the DLS experi- 
ment ([Porter et ami997l : lRoche et aUllOSQl ) in Berkeley. 
For a historical account o f lepton pair production mea- 
surements in general see (ISpechtl . [2007l) . Here we focus 
on the most recent measurements at ultra-relativistic en- 
ergies and the current status of their interpretation. To 
search for non-trivial contributions the data for di-lepton 
measurements are compared to predictions for yields re- 
sulting from the electromagnetic decay of hadrons at 
chemical freeze-out. The hadron production rate is ei- 
ther directly measured or inferr ed from statistical mode l 
calcula tions discussed above ( Braun-Munzinger et all . 
l2004al ). The resulting yields are called 'hadronic cock- 
tail' as they result from the standard known mixture of 
unmodified hadronic resonances. 

Pioneering results on the p roduction o f e~^e~ - 
pairs came from the DL S ([Porter et al. I, ll99<D . 
([Angelis et al. 



HELIOS 



199i 



and CERES 



([Agakichiev et all ll995H200a: lAgakishiev et all . 1 19981 ) 
collaborations: the main and dramatic outcome of 
these experiments was that all central nucleus-nucleus 
collision measurements exhibited a yield that is strongly 
enhanced compared to predictions for cocktail decays 
in the invariant mass range 0.2 < rrie+e- < 1-1 GeV. 
Theoretical an alysis of the excess observ ed in the 
CERES data ([Rapp and Wambach[ . [2000f) indicated 
that the enhancement is due to a strong increase of the 
p-meson width in the hot and dense medium formed in 
the collisio n. The excess dis a ppear s for more peripheral 
collisions ( Agakichiev et all . [20051 ) (which exhibit fea- 
tures more like nucleon-nucleon collisions) but for SPS 
energies the be am energy dependence of the observed 
effect is small ( Adamova et all [2003[ ). A satisfactory 



explanation of the excess observed by DLS at much 
lower energies remained missing. 

Dramatic progress was recently achieved by the NA60 
collaboration which, for a coll ision system of interm ediate 
mass (In-f-In), provided data ( Arnaldi et aL[ . [2b06[ ) in the 
di-muon channel with very good statistics and improved 
mass resolution compared to previous measurements. 
The quality of the data is such that the di-lepton yield 
resulting from final state hadron decays, i.e. the cock- 
tail yield, can be subtracted from the measured di-lepton 
spectra. The resulting subtracted spectrum is compared, 
in Fig.[T21 with predictions that take into account all col- 
lision processes of the p meson with the surrounding par - 
ticles in the fireball ([van Hees and Rappl . l2006l l2007b[ ). 
Note that the data are not acceptance corrected, imply- 
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FIG. 12 The NA60 low mass spectrum after subtrac- 
tion of 'cocktail' decay processes after thermal freeze- 
out compared with th e theoretical predictions from 
refs. ([van Hees and RappL l2006l . l2007bl ). The open data 
points exhibit the size of the correction resultinpj from a de- 
crease in the rj yield by 10% ([Arnaldi et all 12006 



ing that the calculations have to be filtered appropri- 
ately for a m eaningful comparison ([Arnaldi et all . 120061 : 



ISDechtl . [2007h . 

In such many-body calculations, the spectral function 
of the p meson is considerably broadened in the hot and 
dense medium compared to the line shape of the p me- 
son in vacuum. The strong broadening is dominantly 
due to interactions of p mesons with baryons (and anti- 
baryons) in the dense fireball near the phase boundary. 
This is indeed observed in the NA60 data, as is demon- 
strated by the quantitative agreement between data and 
calculations. Note that there is no evidence for a possi- 
ble down ward shift of the p mas s, as had been predicted 
early-on ( Brown and Rhd . Il99lh based on a scaling rela- 
tion ()15p between the p mass and the in-medium quark 
condensate. 

The CERES colla b oratio n has recently pre- 
sented ( Adamova et MI, l2006[ ) their absolutely nor- 
malized data on low-mass e^e~-pair production in 
central Pb-Au collisions at SPS energy, taken with 
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the upgraded CERES apparatus. Again, to explicitly 
display the shape of the in-medium contribution to the 
di-lepton mass spectrum, the cocktail excluding the 
p-meson contribution was subtracted from both the 
data and theoretical calculations. The result is shown 
in Fig. M 
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FIG. 13 The CERES low mass spectrum compared 
with t heoretical predict i ons fr om the 'dropping mass' sce- 
nario (iBrown and Rhd. Il99ll) an d hadronic many-body 
theor y (jvan Hees and Rap d, l2007al lbl: iRapp and Wambachl . 
I2OOOI '). 



Note that the yield of the 'cocktail p' is only about 10 
% of the observed yield near the mass of the vacuum-p 
, demonstrating the extremely strong modification of its 
spectral function in the dense firebal l. Calculations based 
on th e many-body approach of ( van Hees and Rapd . 



l2007al lbl: lRapp and Wambachl . |2000( ) explain this medium 
modification quantitatively, wh ile those based on a 
downward shift of the p mass (|Brown and Rhd . llQQlt 
Ivan Hees and Rappl . l2007al) are at variance with the ob- 
servations^^, in accord with the findings of the NA60 
collaboration. Below di-lepton masses of 200 MeV, the 
CERES data indicate a further strong rise. Such an 
increase towards the 'photon point' (mg+g- = 0) was 
predicted in a consistent treatment of the in-mediu m p- 
meson spectral function ( Rapp and Wambachl . l200dl ) and 
its observation lends further support to the underlying 
theoretical approach. 

Two more experiments have released data on di- 
lepton production in nucleus-nucleus collisions during 
the past year. The HADES collaboration presented 
their first data on C-l-C colli si ons at relatively l ow en - 
ergy ( Agakichiev et aZl . 120071 : lAgakishiev et al. I, [2003), 



substantially corroborating the measurements of the DLS 
collaboration. Currently a significant theoretical effort is 
underway to understand these observations. 

The PHENIX experiment at RHIC has also pre- 
sented firs t results on di- l epton production in Au-Au 
collisions ( Afanasiev et all |2007| ) at very high energy 
{^/snn = 200 GeV). The results are presented in Fig. [HI 
In addition to peaks of the vector mesons, one observes a 
very large enhancement compared to the hadronic cock- 
tail in the di-lepton yield for masses between 200 and 
800 MeV. At present, the size of this enhancement is not 
reproduced within the theoretical approaches described 
above. Future research will tell whether new physics is 
visible here or whether these data can also be described 
within the language of hadronic many-body theories. 



as Au+Au at \[sj^ = 200 GeV 

J/T -» ee 
ee 

cc -> ee (PYTHIA) 



cc -» ee {random correlation) 




See l|Brown et al I l2007t ) for an updated view on the connection 
between the p mass and the chiral condensate. 
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FIG. 14 The PHE NIX low-mass spectrum for Au-Au colli- 
sions at top RHIC (|Afanasiev et all 120071 ') . 

Finally we would like to comment on the connection 
between the medium modification of hadrons and the 
phase diagram. Already 25 years ago Pisarski ( Pisarskil 
Il982[ ) argued that prompt di-lepton production from a 
hot fireball can be a signal for critical behavior through 
changes of the mass and width of the p meson near 
the phase bound a ry. T he issue was further investigated 
in ( Karsch et all Il993( ). As discussed above, it is our 
current understanding that near (or at) the deconfine- 
mcnt phase boundary also chiral symmetry will be re- 
stored. The in-medium electromagnetic response, which 
is dominated by the vector mesons p, uj and 4>, provides 
a direct link to chiral symmetry and its restoration near 
Tc- Restoration of chiral symmetry implies a strong re- 
duction ('melting') of the quark condensate near Tg. Fur- 
thermore, at the phase boundary, the vector- and axial- 
vector correlation functions (see above) corresponding to 
the p meson and its chiral (parity) partner, the ai me- 
son, mu st become identical in the limit of vanishing quark 
masses (jKapusta and Shurvak . 1994; W einberg , 1 9671). 

It is an interesting observation (jRapp and Wambachl . 
I1999I ) that the yield calculated using the hadronic in- 
medium correlation function near the phase boundary co- 
incides remarkably well with that obtained from lowest- 
order qq annihilation in the QGP, where chiral symmetry 
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is restored. Since a strong increase of the p-meson width 
is seen in the present di-lepton data, it thus seems that 
the signal for chiral symmetry restoration in the electro- 
magnetic response of hot and dense matter is a smooth 
'melting' of the p meson i rito a featureless quark- 
antiquark conti nuum (see also ( Gallmeister et all . I2OOII : 
lKamDfeil . [2007l) '). While this is not a rigorous argument 
for chiral symmetry restoration by itself, a much stronger 
case could be made if the modification of the 'chiral' part- 
ner, the Oi could be measured. On general grounds its 
spectral distribution has to become degenerate with that 
of the p meson when chiral symmetry is restored. Hence, 
also the ai meson has to melt smoothly into a quark- 
antiquark continuum. Unfortunately this is hard to check 
experimentally since the dominant electromagnetic decay 
of the oi meson involves, besides a virtual or real pho- 
ton, a pion. The latter suffers strong rescattering and 
absorption in the fireball and hence the early stages of 
the collision are hard to probe. At low temperatures and 
densities, however, one can prove rigorously that chiral 
symmetry restoration manifests itself in a mixing of the 
p- and ai meson through the absorp tion or emission o f a 
pion from the surrounding medium ( Dev all . Il990l ). 



D. Quarkonia-messengers of deconfinement 

In a nucleus-nucleus collision at very high energy heavy 
quarks, charm or beauty, can be produced rather copi- 
ously. For example, the number of charm and anti-charm 
quark pairs in a Pb-Pb collision at LHC energy might 
well reach beyond a hundred. Because of the large mass 
of the charm quarks compared with the typical QCD 
scale {h-QCD ~ 200 McV, mcharm ~ 1-3 GeV) there is 
a separation of time scales between charm quark pro- 
duction and the production of hadrons containing charm 
quarks ( Andronic et aLl . l2008b[ ). The question of medium 
modifications of such hadrons is then more subtle. 

Particles collectively known as 'quarkonia' are bound 
states of charm or beauty quarks and their antiquarks. 
They play an important role as probes for deconfined 
matter inside the hot and de nse fireball. In their sem inal 
1986 paper Matsui and Satz ( Matsui and Satj . ll986t ) ar- 
gued that the bound state made up of charmed quarks 
and anti-quarks, the 3 ftp meson, would be destroyed (or 
prevented from being formed) by the high density of par- 
tons in the QGP. The physics behind this process is sim- 
ilar to Debye screening of the electromagnetic field in 
an electromagnetic plasma through the presence of mov- 
able electric charges. To provide a first estimate, we 
note that the density of partons (quarks and gluons) in 



Whether the melting occurs only close to the phase boundary or 
over an extended range of temperatures and densities below the 
phase boundary is currently an open issue. There are indications, 
however, that already at chemical freeze-out the b aryon density 
is rather low llBraun-Munzinger and Stachell . |2002|) . 



a non-interacting plasma with three massless flavors is 
n = 4.2r^. At a temperature of 500 MeV, this implies 
that n « 70/fm'^. The mean distance between these color 
charges scales like cx 1/T and is about 0.25 fm 

in the ideal gas limit, much less than the spatial extent 
of the J/ip meson. Indeed, taking strong interactions 
among the color charges into account leads to a "Debye 
screening" radius rjy oc l/{gs{T)T) which decreases with 
increasing temperature. Hence the resulting color screen- 
ing may destroy the bound state. The suppressed yields 
of charmonia measured in a high-energy nucleus- nucleus 
collision (compared to their production in the absen ce of 
a QGP) was thus proposed ( Matsui and Satj . ll986[ ) as a 
'smoking gun' signature for the QGP. 

Measurements performed during the last decade at the 
CERN SPS acceler ator indeed pr ovided first evidence for 
such a suppression (lAbreul . [200lh in central collisions be- 
tween heavy nuclei. Little suppression was found in graz- 
ing collisions or collisions between very light nuclei, where 
QGP formatfo n is not expected. The precision data of 
( Abreul . [2OOII ) could be described, however, also by con- 
sidering "normal" absorption of charmonium in the nu- 
clear medium, in conjunction with its possible break-up 
by hadrons produced in the collision (co- movers). Such 
mechanisms could lead to charmoni um suppression even 
in the absence of QGP formation (ICapella et all . l2002t 
iGavin and Vogtl . Il997l : ISpieles et all I1999D and the in- 
terpretation of the SPS data remains inconclusive. 

This situation took an interesting turn in 2000, when 
it was realized that the large number of charm-quark 
pairs produced in a nuclear collisions at RHIC or LHC 
energies leads to new mechanisms for charmonium pro- 
duction, either t hrough statistical production at the 
phase boundary ( Braun-Munzinger and Stachell . l2000l 
I2OOII ). o r through coalescen ce of charm quarks in the 
plasma ( Thews et all . l200ll ). At low energy, the mean 
number of charm-quark pairs produced in a collisions is 
much less than 1, implying that a charmonium state, if 
at all, is always formed from charm quarks of the one 
and only pair produced. On the other hand, the number 
of charm quark pairs at RHIC energies is already much 
larger than 1 (indirect measurements imply a charm- 
quark multiplicity of about 10) and the total number 
of charm quarks in a collision at the LHC is expected 
to reach values larger than 100! Under such conditions 
charm quarks from different pairs can combine to form 
charmonium. Charm-quark recombination works effec- 
tively only if the charm quarks can travel a significant 
distance in the plasma to 'meet' with their prospective 
partner. Under these conditions, charmonium produc- 
tion scales quadratically with the number of charm quark 
pairs. Thus enha ncement, rather than st rong suppres- 
sion, is predicted ( Andronic et all . l2007b[ ) for LHC en- 
ergies. We note that, in the recombination model, it 
is assumed that charmonia are either not formed before 
the QGP or that they are completely destroyed by it 
(complete quenching), so that all charmonium produc- 
tion takes place when the charm quarks hadronizc at the 
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pha se boundary. For a det ailed discussion of this point 
see ([Andronic et all l2008bD . 

The most recent data from the RHIC accelerator pro- 
vide interesting new insight into the connection between 
QGP formation and charmonium production but the 
question how to use charmonia as messengers of decon- 
fincmcnt is far from settled. Here wc briefly describe 
the surprizing aspects of the new PHENIX data and 
argue that they lend first support to the regeneration 
scenario described aboye. The major new insight came 
from a study of the rapidity^^ and ccntrality dependence 
(measured through the number of participating nucle- 
ons in the collision) of the nuclear modification factor 
^AA which has, for the fi rst time, been me asured by 



the PHENIX collaboration (jAdare et aLl . l2007t) in Au-Au 
collisions. This modification factor for J/tJj production is 
defined as 



R 



J/4' 
AA 



^^AuAu/^y 



(17) 



and relates the charmonium yield in nucleus- nucleus colli- 
sions to that expected for a superposition of independent 
nucleon-nucleon collisions. Here, dNj/^,/dy is the rapid- 
ity density of the J/ip yield for AA and pp collisions and 
Ncoii is the number of binary nucleon-nucleon collisions 
for a given centrality class. 



the rapidity dependence of R'^^ at mid-rapidity (cor- 
responding to y=0, i.e. production perpendicular to the 
beam direction). This maximum was entirely unexpected 
as the observed trei id is opposite to t hat expected from 
the melting model (|Matsui and Satj . Il986l : ISatd . |2006[) . 

where i?^{f should attain its smallest value (maximum 
suppression) in regions of phase space with maximum 
energy density, i.e. near mid-rapidity. Likewise, the de- 
struction of charmonia by co-moving hadrons would also 
lead to the largest suppression at mid-rapidity, in conflict 
with PHENIX data. 

On the other hand, the observed maximum of R'^^^' at 
midrapidit y is naturally explained in the recombination 
model of ( Andronic et all l2007a[ ) as being due to en- 



hanced charmonium production at the phase boundary: 
the number of charm quarks is maximal at mid-rapidity 
and this maximum is enhanced even further through re- 
combination. This mechanism provides a good descrip- 
tion of the data, as indicated by the calculated curves 
in Fig. [151 In this sense, the PHENIX measurement 
constitutes first evidence for the statistical production 
of J/ip at chemical freeze-out. Further support for this 
interpretation comes from the observed centrality depen- 
dence of at midrapidity as shown in Fig. [161 We 
reiterate that, if the recombination model is correct, it 
implies complete charmonium quenching in the QGP, as 
discussed above. 
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FIG. 15 Rapidity dependence of the nuclear modification fac- 
tor W l/f for two centrality classes. The data (jAdare et all 
l2007l l (symbols with errors) are compared to calculations 
(lines) . The dashed lines result from uncertainties in the J/tp 
distribution in proton-proton collisions. 
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In Fig. [121 we present the new data from the PHENIX 
collaboration ( Adare et all\206% . The most striking fea- 
ture of these data is the observation of a maximum in 



FIG. 16 Centrality dependence of the nuclear modification 
factor R-^J^ for RHIC and LHC e nergies. The RHIC data 
are from the PHENIX experiment (jAdare et a/.l . l2007l l . The 
calculations are performed w ithin the framework of the sta- 
tistical hadronization model (j Andronic et aLll2007al fbl). 



The 'rapidity' of a particle is defined through its total energy 
E and the longitudinal momentum along the beam axis as: 
y = i In ( E^p" ) ' contrast to a particle's velocity its rapidity 
y is additive under Lorentz transformations. 



One should note that, at present, there are also 
other interpretations of the PHENIX data, in particular 
through cold nuclear matter effects possibly reducing the 
number of gluons and, hence, charm quarks when going 
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away from mid-rapidity. Precision data on J /ip produc- 
tion and its possible hydrodynamic flow will be needed 
to distinguish between these diff erent descriptions. The 
situation is concisely reviewed in ( Granier de Cassagnad . 

Since the number of charm quark pairs is still rather 
moderate at RHIC energies, a strong enhancement of 
J/V' production is not expected in the PHENIX mea- 
surement. The model predictions reproduce very well 
the d ecreasing trend ver sus centrality seen in the RHIC 
data fiAdare et a/.l . l200"7[ ). In contrast, at the much higher 
LHC energy, the charm production cross section is ex- 
pecte d to be about an order of magnitude larger (IVoetl . 
120031 ). As a result, a totally opposite trend as a function 
of centrality is predicted (Fig. [16]) with Raa exceeding 1 
for central collisions. 

For these predictions it is assumed that charm quarks 
are effectively thermalized in the very hot and dense 
QGP, implying that their recombination at the phase 
boundary gives rise to a significant increase in yield near 
mid-rapidity. As shown, the resulting predictions for 
measurements at LHC energies lead to a rather dra- 
matic enhancement rather than suppression in central 
Pb~Pb collisions. If observed, this would be a spectac- 
ular fingerprint of a high-energy quark-gluon plasma, in 
which charm quarks are effectively deconfined. The data 
on charmonium production in Pb-Pb collisions from the 
LHC will be decisive in settling the issue. 



V. PHASES AT HIGH BARYON DENSITY 

Nuclear matter can be compressed in two distinctly dif- 
ferent ways: a rapid squeeze that leads to strong heating, 
as realized in heavy-ion collisions at relativistic energies; 
a slow squeeze which results in cold matter at very high 
baryon density. This type of compression is impossible to 
achieve in the laboratory but is realized in the interior of 
a neutron star, a few seconds after it has been born in a 
supernova explosion. It is conceivable that, in the inner 
core, densities as high as ten tim es that in the middle of 
a hea vy nucleus can be reached ( Lattimer and Prakashl 
I2OO7D . Under such conditions it is expected that the 
closely packed neutrons (with a small admixture of pro- 
tons, electrons, and muons as well as baryons carrying 
strange quarks) dissolve into their constituents and the 
u,d,s-quarks form a degenerate Fermi liquid^^. The com- 
position is determined by charge and color neutrality and 
the requirement of (3 equilibrium, i.e. equilibrium of weak 
interaction processes. 

It has long been known that fermionic systems at low 
temperatures become unstable to the formation and con- 
densation of 'Cooper pairs' if the interaction between 



two fermions is attractive. This situation is expected 
in qu ark-gluon matter above t he de confinement transi- 
tion ijBarroisl . Il977l : iFrautschil . [iQli). Here the Cooper 
instability of the Fermi surface and the formation of di- 
quark pairs is mediated by the attractive interaction in- 
duced by gluon exchange between two quarks of specific 
color, flavor and spin combinations. Since such combina- 
tions carry net color, the new state is called a ' color su- 
perconductor' a. term that was first used in refs. ( Barroisl 
I1977I : IFrautschil . I 19781 ). The presence of color supercon- 
ductivity in the core of neutron stars could lead to in- 
teresting new effects in the long-time evolution of such 
objects such as modifications of the coohng rate through 
neutrino emission, instabilities caused by gravitational 
wave radiation of pulsa rs or glitches in the spin-down 
rate (jAlford et aZ.I . l2007t) . 



A. Color Superconductivity 

Early analysis of the possible pairing patterns in cold 
quark matter and estimates of the resulting gaps A based 
on th e exchange of a single gluon ( Bailin and Lovel Il979l 
Il984[ ) led to values of a few MeV for A. Such low values 
have little influence on the high density EoS. This situa- 
tion changed in the 1990's when it was pointed out that, 
in the physically interesting region of ^ 500 MeV 
which corresponds to about ten times the density in a 
heavy nucleus, perturbative one-giuon exchange is inade- 
quate because of the strong increase in as at such momen- 
tum scales. The resulting non-perturbative effects in the 
quark-gluon coupling were estimated in th e NJL mode l 
and led to gap va l ues o f up to 100 MeV ( Alford et al\ . 
Il998l : iRapp e^"all . Il998l ). Subsequently, it was found 
that the many possible combinations of flavor-color and 
spin degrees of freedom, dictated by the fermionic an- 
tisymmetry of the Cooper pair wave fu n ction , can lead 
to a rich phase stru cture ( Alford et alV l2007t iBuballal . 
I2OO5I: iRischkd . 12004) . For total spin S' = one has 
the possibility of pairing two quark flavors, say up or 
down, leaving the third flavor unpaired (the socalled 
'2SC' phase) or all three quark flavors can participate. 
In this case there is a definite combination of color 
and flavor degrees of freedom called 'color-flavor lock- 
ing' (CFL). Under the conditions of charge and color 
neutrality as well as (3 equilibrium the Fermi energies 
of quarks with given color and flavor quantum numbers 
are in general not equal. The imbalance is partly caused 
by the mass difference — niu^d- For a large mis- 
match, pairing with unequal quantum numbers becomes 
difficult ('stressed superco nductivity') and can even lea d 
to 'gap less' phases ( g 2SC (IShovkovv and Huan"3 . |2003[ ). 
gCFL (jAlford et all \2004 i). Also crystalline phases 



Because of their much larger (Higgs) masses charm, bottom and 
top quarks play no role at the relevant densities. 



similar to the LOFF p h ases ( Fulde and Ferrelll . Il964t 
iLarkin and Ovchinnikovl . Il965l) in conventiona l super- 
conductors are conceivable ( Alford et all l2007t ) . Which 
phase is favored at a given temperature and density is 
determined by the global minimum of the free energy. 
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An example is shown in Fig. 1171 Comparing Fig. [TT] and 
Fig. [T7] it is unlikely that any of the high-density and 
low-temperature phases can be explored in heavy-ion col- 
lisions. 
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FIG. 17 Color superconducting phases at high baryo- 
chemical potential fiq = nt/S and low temperatures as 
predi cted by the NJL m odel in the Hartree approxima- 
tion (|Ruester et 0^.1120051 ). The region of spontaneously bro- 
ken chiral symmetry is denoted by xSB while regions where 
quark matter is in the normal state are indicated by NQ. 
The bold solid lines mark boundaries of first-order transi- 
tions, while the thin lines denote second-order boundaries. 
The dashed lines indicate the boundaries between gapless and 
gapped regions. Several critical points are found. 



Stressed superconductivity can however be studied 
experi mentally in trapped ultracold fermionic atomic 
gases (jGiorgini et al Here an imbalance in chem- 

ical potentials can be achieved by populating two hyper- 
fine states of the atom with a different number of parti- 
cles. At the same time the interaction strength can be 
controlled using Feshbach resonances, to drive the sys- 
tem from weak coupling (BCS regime) through the point 
where diatomic bound states form to the point where 
diatomic molecules undergo Bose-Einstein condensation 
(BEC regime). Thus many of the predicted phases of cold 
quark matter can be 'simulated' in the laboratory with 
interesting future perspectives and cross-fertilization. 

Even though the NJL model is useful in exploring the 
many possibilities of superconducting phases, its quanti- 
tative predictive power is limited by the large sensitiv- 
ity of the results to the model parameters. First prin- 
ciple calculations, on the other hand, are very difficult 
since they require accurate knowledge of the di-quark 
interaction on scales of the Fermi energy ep where as is 
large. Only at very high densities or asymptotically large 
fiq, the coupling becomes small enough to make reliable 
predictions from first principles. In this case one-gluon 
exchange between di-quarks dominates. In the dense 
medium its longitudinal (color-electric) component is De- 
bye screened while the transversal (color-magnetic) com- 
ponents are dynamically screened due to Landau damp- 



ing. This implies that the ratio of the magnetic to electric 
polarization functions goes like ui/\q\ where co is the fre- 
quency and q the three-momentum of the gluon field. In 
the static limit a; — the magnetic components therefore 
remain unscreened. As a consequence, in contrast to the 
usual BCS theory where the pairing gap as a function of 
the coupling constant g varies as A/// ^ exp(— const/g^), 
one has (|Son. . il99a) 



A / 3^^ 
— ~ exp -=— 



(18) 



Such retardation effects for long-ran ge forces are also 
know n in condensed matter physics ( EliashberS . Il960l 
Il96l[ ). The l/(7s-dependence in the exponent of the gap 
function leads to the surprizing phenomenon that the 
pairing gap can take ar bitrarily large values, even thou gh 
the coupling decreases ( Raiagopal and Wilczekl . l200l[ )^^. 
Taking into account the color-flavor-spin degrees of free- 
dom one finds the CFL phase to be the energetically most 
favored pairing state at asymptotically large quark chem- 
ical potentials fFig. [T7|) . 

Even though these ab-initio findings are quite inter- 
esting from a many-body point of view, they are valid 
only for asymptotically large values of Hg, because of 
the logarithmic running of as (Eq. Hence, they 

arc of little relevance for the interior of neutron stars 
where ~ 400 — 600 MeV. One can try to remedy 
this by the inclusion of higher-order corrections in g^. 
Since at such scales gs — 1 it is questionable, however, 
whether such perturbative expansion schemes are justi- 
fied. A more promising approach is to use Schwinger- 
Dyson equations where both the quark and gluon fields 
are treated non-perturbatively with a proper treatment 
of infra -red (small-momentu m behavior) of a^- Recent 
results (|Nickel et all l2006allbl ) indicate that, in the rel- 



evant regime of quark densities in the core of neutron 
stars, pairing gaps of the order of 100 MeV can be ex- 
pected, confirming the earlier findings within NJL model 
studies. 



VI. SUMMARY AND CONCLUSIONS 

In the 30 years since the first discussions about the 
phases of QCD and the corresponding phase diagram 
there has been tremendous progress in our understand- 
ing of strongly interacting matter at extreme conditions. 
Large experimental campaigns have been mounted and 
have amassed a wealth of new data and led to a series of 



Note that gs 



x/IttoT according to Eq. ^ behaves Uke 
'1/ hi IMq if one assumes that the momentum scale Q is gov- 
erned by fiq. Inserting this into Eq. I|18| l it is clear that the 
exponential drops more slowly than l/fJ-q. 

For weak coupling theory to apply in QCD fiq has to bo of the 
order of 10* MeV! {Fig.llj. 
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discoveries. Here we have concentrated on aspects rele- 
vant to the QCD phase diagram. In particular, we have 
discussed that for symmetric matter {fih — 0) the chem- 
ical freeze-out temperature can be determined with an 
uncertainty of better than 10% from measured hadron 
abundances. We have further argued that the observed 
temperature behavior lends strong support to the no- 
tion of a critical temperature Th introduced by Hagcdorn 
and provided arguments that Th coincides with Tc, the 
critical temperature for the quark-hadron transition of 
strongly interacting matter. Thus an important point in 
the phase diagram has been established experimentally. 
We have furthermore summarized the evidence for mass 
changes of hadrons near the phase boundary, with par- 
ticular emphasis on the p meson and laid out arguments 
how these findings arc connected to the restoration of 
chiral symmetry near the phase transition line. Charmo- 
nium production is apparently also strongly influenced 
by the QCD phase transition and we have outlined the 
particular role of this production process for studies of de- 
confinement. Along with the experimental progress also 
came impressive theoretical developments, both concern- 
ing "exact" solutions of QCD on a discrete space-time 
lattice as well as the development of powerful, effective 
models to study the physical processes emerging from the 
experimental observations. 

What may be expected in the future? With the exper- 
imental program at RHIC and in particular the heavy- 
ion program at the CERN LHC the structure of the 
matter above Tc and at vanishing chemical potential can 
be studied quantitatively. In particular, the fireballs 
formed in Pb-Pb collisions at LHC energies will have 
much higher initial temperatures, maybe reaching 1 GeV, 
and live much longer (>10 fm lifetime up to the quark- 
hadron phase transition) than those produced at RHIC. 
Furthermore, hard probes, in particular high transverse- 
momentum jets and heavy quarks, will be abundantly 
produced. From studies in this new environment should 
emerge not only detailed tests of ab-initio QCD predic- 
tions about the phase transition as well as information 
about the bulk properties of the QGP at high tempera- 
ture and its stopping power for high momentum quarks 
but also insight into the nature of the processes that 
lead to confinement. Studies of the phases of strongly 
interacting matter at high densities and moderate tem- 
peratures, on the other hand, are still in their infancy. 
The development of relevant effective theories (including 
the complex reaction dynamics) as well as developments 
of lattice QCD simulations at finite chemical potentials 
arc important milestones in the understanding of quark 
matter at high densities. Further experimental studies 
at lower energy at the RHIC collider as well as with the 



The LHC is now scheduled to start operations in the summer of 
2009, first with protons and afterwards with a pilot run for the 
Pb beam program. 



planned CBM experiment at the FAIR facility at GSI are 
mandatory to make progress in our understanding of the 
QCD phase transition in the high density regime. 

Acknowledgments 

We thank A. Andronic, M. Buballa, K. Rcdlich and A. 
Richtcr for a critical reading of the manuscript. 

References 

Abelev, B. I., et al. (STAR), 2007, Phys. Rev. Lett. 98, 
192301. 

Abreu, M. C. e. a. (NA50), 2001, Phys. Lett. B521, 195. 
Adamova, D., et al. (CERES/NA45), 2003, Phys. Rev. Lett. 
91, 042301. 

Adamova, D., et al, 2006, eprint nud-ex/0611022. 
Adams, J., et al. (STAR), 2005a, Phys. Rev. C72, 014904. 
Adams, J., et al. (STAR), 2005b, Nucl. Phys. A757, 102. 
Adare, A., et al. (PHENIX), 2007, Phys. Rev. Lett. 98, 
232301. 

Adcox, K., et al. (PHENIX), 2005, Nucl. Phys. A757, 184. 
Adler, S. S., et al. (PHENIX), 2006a, Phys. Rev. Lett. 97, 
052301. 

Adler, S. S., et al. (PHENIX), 2006b, Phys. Rev. Lett. 96, 
032301. 

Adler, S. S., et al. (PHENIX), 2007, Phys. Rev. C76, 034904. 
Afanasiev, S., et al. (PHENIX), 2007, eprint arXiv:0706.3034 
[nucl-ex] . 

Agakichiev, G., R. Baur, A. Breskin, R. Chechik, A. Drees, 
C. Jacob, U. Faschingbauer, P. Fischer, Z. Fraenkel, 
C. Fuchs, E. Gatti, P. Glassel, et al, 1995, Phys. Rev. 
Lett. 75, 1272. 

Agakichiev, G., et al (CERES), 2005, Eur. Phys. J. C41, 
475. 

Agakichiev, G., et al (HADES), 2007, Phys. Rev. Lett. 98, 
052302. 

Agakishiev, G., et al (CERES/NA45), 1998, Phys. Lett. 
B422, 405. 

Agakishiev, G., et al (HADES), 2007, eprint arXiv:0711.4281 
[nucl-ex] . 

Akiba, Y. (PHENIX), 2006, Nucl. Phys. A774, 403. 
Alford, M., C. Kouvaris, and K. Rajagopal, 2004, Phys. Rev. 

Lett. 92, 222001. 
Alford, M. G., K. Rajagopal, and F. Wilczek, 1998, Phys. 

Lett. B422, 247. 
Alford, M. G., A. Schmitt, K. Rajagopal, and T. Schafer, 

2007, eprint arXiv:0709.4635 [hep-ph]. 
Allton, C. R., et al, 2003, Phys. Rev. D68, 014507. 
Andronic, A., F. Beutler, P. Braun-Munzinger, K. Redlich, 

and J. Stachel, 2008a, eprint 0804.4132. 
Andronic, A., P. Braun-Munzinger, K. Redlich, and 

J. Stachel, 2007a, Phys. Lett. B652, 259. 
Andronic, A., P. Braun-Munzinger, K. Redlich, and 

J. Stachel, 2007b, Nucl. Phys. A789, 334. 
Andronic, A., P. Braun-Munzinger, K. Redlich, and 

J. Stachel, 2008b, Phys. Lett. B659, 149. 
Andronic, A., P. Braun-Munzinger, and J. Stachel, 2006, 

Nucl. Phys. A772, 167. 
Angelis, A. L. S., et al (HELIOS-3), 1998, Eur. Phys. J. C5, 

63. 



20 



Aoki, Y., G. Endrodi, Z. Fodor, S. D. Katz, and K. K. Szabo, 

2006a, Nature 443, 675. 
Aoki, Y., Z. Fodor, S. D. Katz, and K. K. Szabo, 2006b, Phys. 

Lett. B643, 46. 
Applegate, J. H., C. J. Hogan, and R. J. Scherrer, 1988, ApJ 

329, 572. 

Arnaldi, R., et al. (NA60), 2006, Phys. Rev. Lett. 96, 162302. 
Arsene, I., et al. (BRAHMS), 2005, Nucl. Phys. A757, 1. 
Back, B. B., et al, 2005, Nucl. Phys. A757, 28. 
Baihn, D., and A. Love, 1979, J. Phys. A12, L283. 
Baihn, D., and A. Love, 1984, Phys. Rept. 107, 325. 
Bardeen, J., L. N. Cooper, and J. R. Schrieffer, 1957, Phys. 

Rev. 108, 1175. 
Barrois, B. C, 1977, Nucl. Phys. B129, 390. 
Becattini, F., 1996, Z. Phys. C69, 485. 

Becattini, F., P. Castorina, J. Manninen, and H. Satz, 2008, 

eprint 0805.0964. 
Becattini, F., M. Gazdzicki, A. Keranen, J. Manninen, and 

R. Stock, 2004, Phys. Rev. C69, 024905. 
Becattini, F., and U. W. Heinz, 1997, Z. Phys. C76, 269. 
Becattini, F., J. Manninen, and M. Gazdzicki, 2006, Phys. 

Rev. C73, 044905. 
Bethke, S., 2007, Prog. Part. NucL Phys. 58, 351. 
Blaizot, J.-P., A. Ipp, and A. Rebhan, 2006, Annals Phys. 

321, 2128. 

Boyanovsky, D., H. J. de Vega, and D. J. Schwarz, 2006, Ann. 
Rev. Nucl. Part. Sci. 56, 441. 

Braun-Munzinger, P., K. Redlich, and J. Stachel, 2004a, Par- 
ticle production in heavy ion collisions (World Scientific, 
Singapore), eprint nucl-th/0304013. 

Braun-Munzinger, P., and J. Stachel, 1996a, Nucl. Phys. 
A606, 320. 

Braun-Munzinger, P., and J. Stachel, 1996b, Nucl. Phys. 
A606, 320. 

Braun-Munzinger, P., and J. Stachel, 1998, Nucl. Phys. 
A638, 3. 

Braun-Munzinger, P., and J. Stachel, 2000, Phys. Lett. B490, 
196. 

Braun-Munzinger, P., and J. Stachel, 2001, Nucl. Phys. 
A690, 119. 

Braun-Munzinger, P., and J. Stachel, 2002, Journal ol Physics 

G: Nuclear and Particle Physics 28, 1971. 
Braun-Munzinger, P., and J. Stachel, 2007, Nature 448, 302. 
Braun-Munzinger, P., J. Stachel, J. P. Wessels, and N. Xu, 

1996, Phys. Lett. B365, 1. 
Braun-Munzinger, P., J. Stachel, and C. Wetterich, 2004b, 

Phys. Lett. B596, 61. 
Braun-Munzinger, P., and J. Wambach, 2006, Physik J. 5, 

41. 

Brown, G. E., 1988, Nucl. Phys. A488, 689c. 

Brown, G. E., J. W. Hoh, C.-H. Lee, and M. Rho, 2007, Phys. 

Rept. 439, 161. 
Brown, G. E., and M. Rho, 1991, Phys. Rev. Lett. 66, 2720. 
Buballa, M., 2005, Phys. Rept. 407, 205. 
Cabibbo, N., and G. Parisi, 1975, Phys. Lett. B59, 67. 
Capella, A., A. B. Kaidalov, and D. Sousa, 2002, Phys. Rev. 

C65, 054908. 

Casalderrey-Solana, J., E. V. Shuryak, and D. Teaney, 2005, 

J. Phys. ConL Ser. 27, 22. 
Granier de Cassagnac, R., 2008, eprint 0806.0046. 
Chandrasekhar, S., 1931, Astrophys. J. 74, 81. 
Chaplin, M., 2007, www.lsbu.ac.uk/water/phase.html. 
Cheng, M., et al, 2006, Phys. Rev. D74, 054507. 
Chodos, A., R. L. Jaffe, K. Johnson, C. B. Thorn, and V. F. 



Weisskopf, 1974, Phys. Rev. D9, 3471. 
Cleymans, J., and K. Redlich, 1998, Phys. Rev. Lett. 81, 
5284. 

Collins, J. C, and M. J. Perry, 1975, Phys. Rev. Lett. 34, 
1353. 

Dey, M., V. L. Eletsky, and B. L. loffe, 1990, Phys. Lett. 
B252, 620. 

Dokshitzer, Y. L., and D. E. Kharzeev, 2001, Phys. Lett. 
B519, 199. 

Eliashberg, G. M., 1960, Sov. Phys. JETP 11, 696. 

Eliashberg, G. M., 1961, Sov. Phys. JETP 12, 1000. 

Feinberg, E. L., 1976, Nuovo Cim. A34, 391. 

Fischer, C. S., 2006, Journal of Physics G: Nuclear and Par- 
ticle Physics 32, R253. 

Fisher, M. E., 1967, Physics 3, 255. 

Fodor, Z., and S. D. Katz, 2002, Phys. Lett. B534, 87. 

Fodor, Z., and S. D. Katz, 2004, JHEP 04, 050. 

de Forcrand, P., and O. Philipsen, 2002, Nucl. Phys. B642, 
290. 

Fortov, V. E., R. I. Ilkaev, V. A. Arinin, V. V. Burtzev, 
V. A. Golubev, I. L. losilevskiy, V. V. Khrustalev, A. L. 
Mikhailov, M. A. Mochalov, V. Y. Ternovoi, and M. V. 
Zhernokletov, 2007, Physical Review Letters 99, 185001 
(pages 4). 

Frautschi, S. C, 1971, Phys. Rev. D3, 2821. 

Frautschi, S. C, 1978, in Proc. of the Workshop on Hadronic 

Matter at Extreme energy density, edited by N. Cabbibo 

(Erice, Italy). 

Fritzsch, H., M. Gell-Mann, and H. Leutwyler, 1973, Phys. 

Lett. B47, 365. 
Fromerth, M. J., and J. Rafelski, 2002, eprint astro- 

ph/0211346. 

Fulde, P., and R. A. FerreU, 1964, Phys. Rev. 135, A550. 
Gallmeister, K., B. Kampfer, O. Pavlenko, and C. Gale, 2001, 

Nucl. Phys. A688, 939. 
Gavin, S., and R. Vogt, 1997, Phys. Rev. Lett. 78, 1006. 
Gebhard, F., 1997, The Mott Metal- Insulator Transition 

(Springer, Berlin). 
Giorgini, S., L. P. Pitaevskii, and S. Stringari, 2007, Theory 

of ultracold fermi gases, eprint arXiv:0706.3360. 
Goldstone, J., 1961, Nuovo Cim. 19, 154. 
Greenberg, O. W., 1964, Phys. Rev. Lett. 13, 598. 
Gross, D. J., and F. Wilczek, 1973, Phys. Rev. Lett. 30, 1343. 
Gyulassy, M., and L. McLerran, 2005, Nucl. Phys. A750, 30. 
Hagedorn, R., 1965, Nuovo Cim. Suppl. 3, 147. 
Heckmann, K., 2007, private communication, 
van Hees, H., and R. Rapp, 2006, Phys. Rev. Lett. 97, 102301. 
van Hees, H., and R. Rapp, 2007a, private communication, 
van Hees, H., and R. Rapp, 2007b, eprint arXiv:0711.3444 

[hep-ph]. 

Heinz, U., and G. Kestin, 2006, PoS CPOD2006, 038. 
Heinz, U. W., 1998, Nucl. Phys. A638, 357c. 
Huang, K., and S. Weinberg, 1970, Phys. Rev. Lett. 25, 895. 
Huovinen, P., P. F. Kolb, U. W. Heinz, P. V. Ruuskanen, and 

S. A. Voloshin, 2001, Phys. Lett. B503, 58. 
Itoh, N., 1970, Prog. Theor. Phys. 44, 291. 
Kampfer, B., 2007, private communication. 
Kampfer, B., M. Bluhm, H. Schade, R. Schulze, and D. Seipt, 

2007, eprint arXiv:0708.3322 [hep-ph]. 
Kapusta, J. I., and E. V. Shuryak, 1994, Phys. Rev. D49, 

4694. 

Karsch, F., 2002, Lect. Notes Phys. 583, 209. 
Karsch, F., K. Redlich, and L. Turko, 1993, Z. Phys. C60, 
519. 



21 



Kolb, P., and U. Heinz, 2004, in Quark-Gluon PlasmaS, edited 
by R. Hwa and X. Wang (World Scientific, Singapore), pp. 
634-714. 

Larkin, A. I., and Y. N. Ovchinnikov, 1965, Sov. Phys. JETP 
20. 

Lattimer, J. M., and M. Prakash, 2007, Phys. Rept. 442, 109. 
Matsui, T., and H. Satz, 1986, Phys. Lett. B178, 416. 
Mott, N. F., 1968, Rev. Mod. Phys. 40, 677. 
Mrowczynski, S., and M. H. Thoma, 2007, Annual Review of 

Nuclear and Particle Science 57, 61. 
Nambu, Y., 1960, Phys. Rev. Lett. 4, 380. 
Nambu, Y., and G. Jona-Lasinio, 1961a, Phys. Rev. 122, 345. 
Nambu, Y., and G. Jona-Lasinio, 1961b, Phys. Rev. 124, 246. 
Nickel, D., R. Alkofer, and J. Wambach, 2006a, Phys. Rev. 

D74, 114015. 

Nickel, D., J. Wambach, and R. Alkofer, 2006b, Phys. Rev. 
D73, 114028. 

O'Hara, K. M., S. L. Hemmer, M. E. Gehm, S. R. Granade, 
and J. E. Thomas, 2002, Science 298, 2179. 

Philipsen, O., 2006, PoS LAT2005, 016. 

Philipsen, O., 2007, eprint arXiv:0710.1217 [hep-ph]. 

Pisarski, R. D., 1982, Phys. Lett. BllO, 155. 

Politzer, H. D., 1973, Phys. Rev. Lett. 30, 1346. 

Porter, R. J., et al. (DLS), 1997, Phys. Rev. Lett. 79, 1229. 

Rafelski, J., and J. Letessier, 2008, J. Phys. G35, 044042. 

Rajagopal, K., and F. Wilczek, 2001, At the frontier of Parti- 
cle Physics/ Handbook of QC'D (Word Scientific), pp. 2061- 
2151. 

Rapp, R., T. Schafer, E. V. Shuryak, and M. Velkovsky, 1998, 

Phys. Rev. Lett. 81, 53. 
Rapp, R., and J. Wambach, 1999, Eur. Phys. J. A6, 415. 
Rapp, R., and J. Wambach, 2000, Adv. Nucl. Phys. 25, 1. 
Proc. Quark-Matter 2005 Conference, 2006, Nucl. Phys. 

A774. 

Proc. Quark-Matter 2006 Conference, 2007, J. Physics G: Nu- 



clear and Particle Physics 34. 

Proc. Quark-Matter 2008 Conference, 2008, J. Physics G: Nu- 
clear and Particle Physics 35. 

Proceedings of Science, 2006, CPOD2006. 

Rischke, D. H., 2004, Prog. Part. Nucl. Phys. 52, 197. 

Roche, G., et al. (DLS), 1989, Phys. Lett. B226, 228. 

Ruester, S. B., V. Werth, M. Buballa, I. A. Shovkovy, and 
D. H. Rischke, 2005, Phys. Rev. D72, 034004. 

Satz, H., 2006, Journal of Physics G: Nuclear and Particle 
Physics 32, R25. 

Schwarz, D. J., 1998, Mod. Phys. Lett. A13, 2771. 

Shapiro, S. L., and S. A. Teukolsky, 1983, Black Holes, White 
Dwarfs and Neutron Stars (John Wiley & Sons, New York). 

Shovkovy, I., and M. Huang, 2003, Phys. Lett. B564, 205. 

Shuryak, E. V., 1978a, Phys. Lett. B78, 150. 

Shuryak, E. V., 1978b, Phys. Lett. B78, 150. 

Son, D. T., 1999, Phys. Rev. D59, 094019. 

Specht, H. J., 2007, eprint arXiv:0710.5433 [nucl-ex]. 

Spieles, C., et al, 1999, Phys. Rev. C60, 054901. 

Stephanov, M. A., 2004, Prog. Theor. Phys. Suppl. 153, 139. 

Stock, R., 1999, Phys. Lett. B456, 277. 

Stoecker, H., 2005, Nucl. Phys. A750, 121. 

Teaney, D., J. Lauret, and E. V. Shuryak, 2002, Nucl. Phys. 
A698, 479. 

Thews, R. L., M. Schroedter, and J. Rafelski, 2001, Phys. 

Rev. C63, 054905. 
Thomas, D., D. N. Schramm, K. A. Olive, and B. D. Fields, 

1993, Astrophys. J. 406, 569. 
Vitev, I., 2006, J. Phys. Conf. Ser. 50, 119. 
Vogt, R. (Hard Probe), 2003, Int. J. Mod. Phys. E12, 211. 
Weinberg, S., 1967, Phys. Rev. Lett. 18, 507. 
Yang, C.-N., and R. L. Mills, 1954, Phys. Rev. 96, 191. 
Yao, W.-M., et al. (Particle Data Group), 2006, Journal of 

Physics G: Nuclear and Particle Physics 33, 1. 



